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These lecture notes give a short review of methods such as the matrix ansatz, the additivity 
principle or the macroscopic fluctuation theory, developed recently in the theory of non-equilibrium 
phenomena. They show how these methods allow to calculate the fluctuations and large deviations 
of the density and of the current in non-equilibrium steady states of systems like exclusion processes. 
The properties of these fluctuations and large deviation functions in non-equilibrium steady states 
(for example non-Gaussian fluctuations of density or non-convexity of the large deviation function 
which generalizes the notion of free energy) are compared with those of systems at equilibrium. 



I. INTRODUCTION 



The goal of these lectures, delivered at the Newton Institute in Cambridge for the workshop " Non-EquiUbrium 
Dynamics of Interacting Particle Systems" in March- April 2006, is to try to introduce some methods used to study 
non-equilibrium steady states for systems with stochastic dynamics and to review some results obtained recently on 
the fluctuations and the large deviations of the density and of the current for such systems. 

Let us start with a few examples of non-equilibrium steady states: 

1. A system in contact with two heat baths at temperatures Ta and Tfj. 
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Figure 1: A system in contact with two heat baths at temperatures Ta and Ti,. 

At equilibrium, i.e. when the two heat baths are at the same temperature (Tq = = T), the probability P{C) 
of finding the system in a certain microscopic configuration C is given by the usual Boltzmann-Gibbs weight 

-Pequilibrmm(C) = -^"^ SXp 

where E{C) is the internal energy of the system in configuration C. Then the task of equilibrium statistical 
mechanics is to derive macroscopic properties (equations of states, phase diagrams, fluctuations,...) from ([1]) as 
a starting point. A very simplifying aspect of ([T]) is that it depends neither on the precise nature of the couplings 
with the heat baths (at least when these couplings are weak) nor on the details of the dynamics. 

When the two temperatures Ta and T;, are different, the system reaches in the long time limit a non-equilibrium 
steady state (T], H, |3, [l] ) but there does not exist d, @ an expression which generalizes ^ for the steady state 
weights P{C) of the microscopic configurations. 
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In fact for a non-equilibrium system, the steady state measure P{C) depends in general on the dynamics of the 
system and on its couplings with the heat baths. 

Beyond trying to know the steady state measure P(C), which can be done only for very few examples 0, H, 
B [3j [O, Ii3.113| , one might wish to determine a number of properties of non-equilibrium steady states like 
the temperature or energy profiles [3, [IE] , the average flow of energy through the system [IB, [13, [H, [I3| , the 
probability distribution of this energy flow, the fluctuations of the internal energy or of the density. 

2. A system in contact with two reservoirs of particles at densities pa and pb- 

Another non-equilibrium steady state situation one can consider is that of a system exchanging particles with 
two reservoirs [2(l| at densities pa and pb- When pa ^ Pb (and in absence of external field) there is a flow of 
particles through the system. One can then ask the same questions as for the previous case: for example what 
is the average current of particles between the two reservoirs, what is the density profile through the system, 
what are the fluctuations or the large deviations of this current or of the density. 
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Figure 2; A system in contact with two reservoirs at densities pa and pb- 
3. The symmetric simple exclusion process (SSEP) 
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Figure 3: The symmetric simple exclusion process. 

The SSEP [m, [13, [13, [13] is one of the simplest models of a system maintained out of equilibrium by contact 
with two reservoirs at densities pa and pb- The model is defined as a one-dimensional lattice of L sites with 
open boundaries, each site being either occupied by a single particle or empty. During every infinitesimal time 
interval dt, each particle has a probability dt of jumping to its left neighboring site if this site is empty, and 
a probability dt of jumping to its right neighboring site if this right neighboring site is empty. At the two 
boundaries the dynamics is modified to mimic the coupling with reservoirs of particles: at the left boundary, 
during each time interval dt, a particle is injected on site 1 with probability adt (if this site is empty) and a 
particle is removed from site 1 with probability ^dt (if this site is occupied). Similarly on site L, particles are 
injected at rate 5 and removed at rate (3. 

We will see ((02]) below and [1^, [1^, H^]) that these choices of the rates a, 7, (3, 6 correspond to the left boundary 
being connected to a reservoir at density pa and the right boundary to a reservoir at density pb with pa and pb 
given by 

a S 
P« = — ] — ; Pb ITT^ ■ (2) 



a + 7 



f3 + S 



One can also think of the SSEP as a simple model of heat transport, if one interprets the particles as quanta of 
energy. Then if each particle carries an energy e, the SSEP becomes the model of a system in contact with two 
heat baths at temperatures Ta and Tb given by (see next section) 
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Figure 4: The asymmetric simple exclusion process. 



4. Driven diffusive systems 

One can add to the systems described above an electric or a gravity field which tends to push the particles in a 
preferred direction. 

For example adding a field to the SSEP means that the hopping rates to the left become q (the hopping rates 
to the right being still 1). The model becomes then the ASEP (the asymmetric simple exclusion process) 
9ll28l . I29L I30L I31I which appears in many contexts , such as hopping conductivity [34| , models of trafiic 

35| , growth [36[ or polymer dynamics [37[ . In presence of this external field, the system reaches a non-equilibrium 
steady state even for a ring geometry, without need of a reservoir. 

The large scale of the ASEP differs noticeably from the SSEP. For example in the ASEP on the infinite line, one 
can observe shock waves whereas the SSEP is purely diffusive. In fact on large scales the ASEP is described (s^ 
by the Kardar Parisi Zhang equation [3HI while the SSEP is in the universalily class of the Edwards Wilkinson 
equation [39l.[40l|. 



The outline of these lectures is as follows: 

In section In] it is recalled how detailed balance should be modified to describe systems in contact with several heat 
baths at unequal temperatures or several reservoirs at different densities. 

In section IIIII the large deviation functional of the density is introduced and there is a comparison between its 
properties in equilibrium and in non-equilibrium steady states. 

In section llVi the connexion between the non-locality of the large deviation functional of the density and the 
presence of long range correlations is discussed. 

In section |V] it is shown how to write the evolution equations of the profile and of the correlation functions for the 
symmetric simple exclusion process. 

Section IVII describes the matrix ansatz [lol | which gives an exact expression of the weights in the non-equilibrium 
steady state of the symmetric exclusion process. 

Using an additivity relation IVIII established as a consequence of the matrix ansatz, the large deviation functional 
[lEHi] of the density for the SSEP is calculated in lVlIIl 

The macroscopic fluctuation theory of Bertini, De Sole, Gabrielli, Jona-Lasinio and Landim [4l|, IH, [3] is 
recalled in section IIXI which shows how the calculation of large deviation functional of the density can be formulated 
as an optimisation problem. 

The definition of the large deviation function of the current and the fluctuation theorem [4^, |4^, [13, |H, |4^ are 
recalled in section |X] from which the fluctuation-dissipation theorem for energy or particle currents can be recovered 
(section IXI[) . 

A perturbative approach ,50] to calculate the large deviation function of the current for the SSEP is sketched in 

Kni 

The additivity principle, which predicts the cumulants and the large deviation function of the current, is presented 
in section IXIIII 

The last four sections are devoted to the ASEP: the matrix ansatz for the ASEP is recalled in section IXIVI It is 
shown in section IXVl how to obtain the phase diagram of the TASEP from the matrix ansatz. An additivity relation 
from which one can compute the large deviation funtional of the density [s^. [ssl is established in section IXVII Latsly 
in section rXVIII it is shown that the fluctuations of density are non-Gaussian [54] in the maximal current phase of the 
TASEP. 
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II. HOW TO GENERALIZE DETAILED BALANCE TO NON-EQUILIBRIUM SYSTEMS 



As in non-equilibrium systems, the steady state measure P{C) depends on the coupHngs to the heat baths and on 
the dynamics of the system, each model of a non-equilibrium has to incorporate a description of these couplings and 
of the dynamics (various ways of modeling the effect of heat baths or of reservoirs are described in see 0,111). It 

IS 

often theoretically simpler to represent the effect of the heat baths (or of the reservoirs of particles) by some stochastic 
terms such as Langevin forces corresponding to the temperatures of the heat baths. In practise the dynamics becomes 
a Markov process. 

For a system with stochastic dynamics given by a Markov process (such as the SSEP or mechanical systems with 
heat baths represented by Langevin forces) the evolution is specified by a transition matrix W{C' , C) which represents 
the rate at which the system jumps from a configuration C to a configuration C (i.e. the probability that the system 
jumps from C to C" during an infinitesimal time interval dt is given by W{C' , C)dt). For simplicity, we will limit the 
discussion to the case where the total number of accessible configurations is finite. The probability Pt{C) of finding 
the system in configuration C at time t evolves therefore according to the Master equation 



dPtjC) 
dt 



Y,WiC,C')Pt{C')-W{C',C)PtiC) . (4) 



One can then wonder what should be assumed on the transition matrix W{C',C) to describe a system in contact 
with one or several heat baths (as for example in figure [T]). 

At equilibrium, (i.e. when the system is in contact with a single heat bath at temperature T) one usually 
requires that the transition matrix satisfies detailed balance 

W{C', C) = W{C, C) e-^ . (5) 

This ensures the time reversal symmetry of the microscopic dynamics: at equilibrium (i.e. if the initial condition is 
chosen according to ([1])), the probability of observing any given history of the system {C^, < s < <} is equal to the 
probability of observing the reversed history 

Pro({C„ 0<s<t}) = Pro({Ct_„ < s < t}) . (6) 

Therefore if e is the energy transferred from the heat bath at temperature T to the system, and We{C' , C)dt is the 
probability that the system jumps during dt from C to C by receiving an energy e from the heat bath, one can rewrite 
the detailed balance condition ([5]) as 

W,{C', C) = e^^ W-,{C, C) . (7) 

If detailed balance gives a good description of the coupling with a single heat bath at temperature T, the straight- 
forward generalization of ([7]) for a system coupled to two heat baths at unequal temperatures like in figure [1] is 

M 

W,„^,,(C',C) = e-^-* W-e^,-,,{C,C') (8) 

where ea,£b are the energies transferred from the heat baths at temperatures Ta^Tf, to the system when the system 
jumps from configuration C to configuration C". By comparing with ([7]), this simply means that the exchanges of 
energy with the heat bath at temperature Ta tend to equilibrate the system at temperature Ta and the exchanges 
with the heat bath at temperature tend to equilibrate the system at temperature Tf,. 

For a system in contact with two reservoirs of particles at fugacities Za and Zj,, as in figure[21 the generalized detailed 
balance ^ becomes 

W,^^,,iC',C) = zl^zf W^_,„,_,,(C,C') (9) 

where qa and are the numbers of particles transferred from the two reservoirs to the system when the system 
jumps from configuration C to configuration C . 



From the definition of the dynamics of the SSEP, it is easy to check that it satisfies © with 



One can also check from ([3]) that if one interprets the particles as quanta of energy, ([8]) is satisfied. 



One way of justifying ([8]) is to consider the composite system made up of the system we want to study and of the 
two reservoirs. This composite system is isolated and therefore its total energy £ 



£ = E{C) + Ea + Eb 



(11) 



is conserved by the dynamics. In ([TT|) E{C) is the energy of the system we want to study and Ea,Ei,. are the energies 
of the two reservoirs (for simplicity we assume that the energy of the coupling between the reservoirs and the system 
is small). Whenever there is an evolution step in the dyanmics, the system jumps from the microcopic configuration 
C to the configuration C and the energies of the reservoirs jump from Ea, Eb to E'^, E'^. For the composite system to 
be able to reach the microcanonical distribution and for microcanonical detailed balance to hold one needs that the 
transition rates satisfy 



Sa(Ea) + Su(Eu) S„ (E' ) + Sfc (E' ) 

e ^ Pro({C7, Ea,Eb} ^ {C, E'„ E'b}) = e ^ Pro({C', E'^} ^ {C, Ea,Eb}) 



(12) 



where Sa{Ea) and Sb{Eb) are the entropies of the two reservoirs at energies Ea and Eb- Then if the heat baths are 
large enough, one has 



S{Ea) - S{E',) 



En — E' 



S{Eb) - S{E'b) 



Eb 



E'b 



Ta ' ' Tb 

where Ta and Tb are the (microcanonical) temperatures of the two heat baths and p2p reduces to 



(13) 



Remark: The quantity —-gf in ([8]) is the entropy produced in the reservoirs. In fact, in the theory 

of non-equilibrium phenomena, one can associate to an arbitrary Markov process, defined by transition rates 
W{C' ,C), an entropy production [H, El, [33, Si, [sS, [sll (in the surrounding heat baths) given by 



AS'(C -> C) = log 



W{C',C) 
W{C, C) 



and ([8]) appears as one particular case of this general definition. 



III. FREE ENERGY AND THE LARGE DEVIATION FUNCTION 



At equilibrium the free energy is defined by 

F = -kTlogZ = -kT log 



^exp 

. c 



In this section we are going to see that the knowledge of the free energy gives also the distribution of the fluctuations 
and the large deviation function of the density. This will enable us to extend the notion of free energy to non- 
equilibrium systems by considering the large deviation functional j2ll . [sol [60| of the density. 




Figure 5; For a system of A'' particles in total volume V , the probability Pviji) of having n particles in a large subvolume v is 
given by p4p . 

If one considers a box of volume V containing N particles as in figure [51 the probability Pv{n) of finding n particles 
in a subvolume v located near a position r has the following large v dependence 



Pv{n) ^ exp —V ap ^— ^ 



(14) 
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where af{p) is a large deviation function. Figure [6] shows a typical shape of afr{p) for an homogeneous system (i.e. 
not at a coexistence between different phases) with a single minimum at p = p* where ap{p) vanishes. 




Figure 6: A typical shape of the large deviation function ap{p). The most likely density p* is the value where ap(p) vanishes. 



One can also define the large deviation functional J- for an arbitrary density profile. If one divides a system of 
linear size L into n boxes of linear size I (in dimension d, one has of course n = L'^/l'^ such boxes), one can try to 
determine the probability of finding a certain density profile {pi, P2, ■■•Pn}, i-e. the probability of seeing I'^pi particles 



in the first box, l'^p2 particles in the second box, 
dependence of this probability 



I Pn in the nth box. For large L one expects the following L 



n 



Figure 7: In (|15|l one specifies the densities pi in each box 



Pro(pi, ...p„) exp[-L''j^(pi,p2, ■■■Pn)] (15) 
where J-" is a large deviation function which generalizes ap{p) defined in ([H]) . If one introduces a reduced coordinate 

X 

f = Lx (16) 

and if one takes the limit L oo, I ^ oo with Z <C i so that the number n of boxes becomes large, this becomes a 
functional !F{p{x)) for an arbitrary density profile p{x) 

Pro(p(f)) ^ eM~L''T{p{x))] . (17) 

Clearly the large deviation function af;{p) or the large deviation functional J-{p{x)) can be defined for equilibrium 
systems as well as for non-equilibrium systems. 



6 



For equilibrium systems, one can show that ap{p) is closely related to the free energy: if the volume v is suffi- 
ciently large, for short ranged interactions and in absence of external potential, the large deviation function ap{p) is 
independent of r and its expression is given by 



ap{p) = a{p) 



f{p)-f{p*)-ip-p*)f'{p*) 
kT 



(18) 



where f{p) is the free energy per unit volume at density p and p* = ^. This can be seen by noticing that if 
much larger than the range of the interactions and if v <^ V one has 



Zy{n) Zv-v{N ~ n) 
MAO 



0{v — ) (19) 

where Zv{N) is the partition function of N particles in a volume V and the term exp 0{v~3~) represents the 

interactions between all pairs of particles, one of which is the volume v and the other one in V — v. Then taking the 
log of (|19p and using the fact that the free energy /(p) per unit volume is defined by 



lim 



log ZviVp) _ f{p) 



V 



kT 



(20) 



one gets p8)) . The functional T can also be expressed in terms of f{p): if one considers Vp* particles in a volume 
V = L'^, one can generalize for systems with short range interactions and no external potential 



„ . ^ Zy{vpi)...Zy{vpn) 

Pro pi, ...pn) = exp 

Zv{Vr) 

where v = l'^. Comparing with (jlSp. in the limit L ^ oo, I ^ oo, keeping n fixed gives 

71 

T{p,,p,,...p.^)^—-J2[f{p,)-f{p*)]. 



(21) 



(22) 



In the limit of an infinite number of boxes, this becomes 



T{pix)) 



1 

kf 



dx [fip{x)) - f{p*)] 



(23) 



Thus for a system at equilibrium, the large deviation functional J- is fully determined by the knowledge of the free 
energy f{p) per unit volume. In (I23p . we see that 

• The functional J-' is a, local functional of p{x). 

• It is also a convex functional of the profile p(x) (as the free energy f{p) is a convex function of the density p, 
i.e. f{api + (1 - a)p2) < a/(pi) + (1 - a)/(p2) for < a < 1) . 

• When /(p) can be expanded around p* (i.e. at densities where the free energy /(p) is not singular) one obtains 
also from that the fluctuations of the density profile are Gaussian. In fact if one expands near p* 
and one replaces it into p4p one gets that the distribution of the number n of particles in the subvolume v is 
Gaussian (if v is large enough) 



Pv (n) ~ exp 



exp 



I"{p*) 
2vkT 



[n — vp*)"^ 



and its variance, as predicted by Smoluchowki and Einstein, is given by 



f"{p* 



V kT k{p*) 



where the compressibility k(p) is defined by 



k(p) 



1 dp 
pdp 



(24) 



(25) 



(26) 



(and the pressure p is given as usual by p = — ^^/(y) = P* f'{p*) — f{p*))- 

Note that at a phase transition, /(p) is singular and the fluctuations of density are in general non-Gaussian. 
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• One also knows (by the Landau argument) that, with short range interactions, there is no phase transition if 
the dimension of space is one dimension. 

In contrast to equilibrium systems, one can observe in non-equilibrium steady states of systems such as the ones 
described in figures [1] and [2] 

• The large deviation functional T may be non local. For example in the case of the SSEP, we will see in section 
IVIIII that the functional is given by for pa — Pb small by (sec ([75)1 below). 



^({p(x)}) = / dx 

{Pa - Pbf 



(27) 



[Pa(l - PaW Jo 

where p*{x) is the most likely profile 



dx / dyx{l-y){p{x)^ p*{x)){p{y)- p*{y)) +0{pa- Pb)^ 



p*{x) = (1 - x)pa + Xpb 



(28) 



• For the ASEP, there is a range of parameters where the functional is non convex (see [52, [5J| and section 
KVll below). 

• There are also cases, where in the maximal current phase, the density fluctuations are non Gaussian (see (s^ l 
and section IXVIII below) . 

• In non-equilibriurn systems nothing prevents the existence of phase transitions in one dimension [lol . [ill [sU . 
M, [62„ M, M [M, [63, [68,, M, Jfl, Ji]. 



IV. NON LOCALITY OF THE LARGE DEVIATION FUNCTIONAL OF THE DENSITY AND LONG 

RANGE CORRELATIONS 



A feature characteristic of non-equilibrium systems is the presence of weak long range correlations |72l . l73l . 174 , l75l . 
[tgI . ■ For example for the SSEP, we will see [T^ in next section (|44l45p that for large L the correlation function of 
the density is given for < a; < y < 1 



{p{x)p{y)), = -^-P^^-P^x{l~y) 



(29) 



The presence of these long-range correlations is directly related to the non-locality of the large deviation functional 
J- . This can be seen by introducing the generating function Q{{a{x)}) of the density defined by 



exp [Lg{{a{x)})] 



cxp 



L I a{x)p{x)dx 



(30) 



where (.) is an average over the profile p{x) in the steady state. As the probability of this profile is given by pT|) 
the average in (j30[) is dominated, for large L, by an optimal profile, which depends on a{x), and G is the Legendre 
transform of J- 



g{{a{x)}) = max 



a{x)p{x)dx — T{{p{x)}) 



(31) 



It is clear from ([5T|) that if the large deviation is local (as in then the generating function Q is also local. Now 

by taking derivatives with respect to a{x) one gets that the average profile and the correlation functions are given by 



P*{X) EE {p{x)) 



sg 



Sa{x) 



(32) 



i(a;)=0 



{p{x)piy)),^ {pix)piy)) ~ {pix)){p{y)) 



L 6a{x) 6a{y) 



(33) 



a{x)=0 
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This shows that the non-loacahty of Q is directly related to the existence of long range correlations. 

One can understand the L dependence in (|33p by assuming that the non-local functional Q is can be expanded as 



1 rl 1-1 1-1 

2 



Q{a{x)) — dx A{x) a{x) + dx B{x) a{x) + dx dy C{x,y) a(x) aijj) + 

Jo Jo Jo Jx 

If one conies back to a discrete system of L sites with a number rii of particles on site i, one has 



Lg(a{x)) ~ log 



exp 



(34) 



(35) 



By expanding in powers of the one has 



2—1 i—1 i<j 

Clearly one has 

and comparing and ([55]) one sees that 



a. 



1 



-c 



which leads to (j33[) . A similar reasoning would show that 

1 5''g 



{p{xi)p{x2)-.p{xk))c = 



L'^ 1 Sa{xi)...Sa{xk) 



(36) 



(37) 



(38) 



(39) 



y(x)=0 



This I/L'^^^ dependence of the k point function can indeed be proved in the SSEP [t^]- We see that all the correlation 
functions can in principle be obtained by expanding, when this expansion is meaningful (see [s^ . Issj ] for counter- 
examples), the large deviation function Q in powers of a{x). 

V. THE SYMMETRIC SIMPLE EXCLUSION MODEL 

For the SSEP, the calculation of the average profile or of the correlation functions can be done directly from the 
definition of the model. If = or 1 is a binary variable indicating whether site i is occupied or empty, one can 
write the time evolution of the average occupation (r^) 



d{ri) 
dt 

d{n) 

dt 

d{TL) 

dt 



=Q! - (a -I- 7 l)(ri) + {T2) 

= (r,_i>-2(T,) + (r,+i> for 2<i<L~l 

= {tl-i)^{1 + P + 5){tl)+6 



(40) 



The steady state density profile (obtained by writing that ^^^^ — 0) is [l^ 

, , Pa{L + 1^ - i) + Pb{i - I + 



a+7 f3+d 



(41) 



with Pa and pb defined as in One can notice that for large L, if one introduces a macroscopic coordinate i = Lx, 
this becomes 



(Ti) = p*ix) = (1 - x)pa + Xpb 



(42) 
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and one recovers p8)) . For large L one can also remark that (ri) — > pa and (tl) ^ Pb indicating that pa and 
defined by ([2]) represent the densities of the left and right reservoirs. One can in fact show [2^, [2^ (2^ that the rates 
a, 7, /3, 5 do correspond to the left and right boundaries being connected respectively to reservoirs at densities pa and 
Pb- 

The average current in the steady state is given by 

(J) = (r.(l - r.+i) - r,+i(l - r,)) = {n - n+i) = f° ~ ^ r (43) 

^+ — + - 

This shows that for large L, the current (J) ~ PslzBl jg proportional to the gradient of the density (with a coefficient 
of proportionality which is here simply 1) and therefore follows Fick's law. 

One can write down the equations which generalize (j40p and govern the time evolution of the two-point function 
or higher correlations. For example one finds [tI, [T^I in the steady state for 1 < i < j < L 



(7^i^ + ^ + ^-i)'(^ + ^ + ^-2) 



For large L, if one introduces macroscopic coordinates i — Lx and j = Ly, this becomes for x < y 

{TL.TLy)c^^^^^^[Pa-pb? (45) 

which is the expression ([29|) . 

One could believe that these weak, but long range, correlations play no role in the large L limit. However if one 
considers macroscopic quantities such as the total number N of particles in the system, one can see that these two 
point correlations give a leading contribution to the variance of N 



{N^) - {Nf = ^[(r.) - {nf] + 2^(r,T,), ^ l\ j dxp*{x){l - p*{x)) - 2(p, ~ p,)' f dx f 



dy x{l - y) 



(46) 



For the SSEP, one can write down the steady state equations satisfied by higher correlation functions to get for 
example for x < y < z 

— 2u)(l — z) , ,0 ,,„, 

{TLxTLyTL.)c = -2^ '-{pa - Pbf (47) 

but solving these equations become quickly quickly too complicated. We will see in next section that the matrix 
ansatz gives an algebraic procedure to calculate all these correlation functions [77| . 



VI. THE MATRIX ANSATZ FOR THE SYMMETRIC EXCLUSION PROCESS 

The matrix ansatz is an approach inspired by the construction of exact eigenstates in quantum spin chains [tsI. ItqI . 
1131 . It gives an algebraic way of calculating exactly the weights of all the configurations in the steady state. In [id] it 
was shown that the probability of a microscopic configuration {Ti,r2, ...tl} can be written as the matrix element of 
a product of L matrices 



PrOll Ti , Tq, ...Tr 1) = — , r , , (48) 

^' ^' {W\{D + E)L\V) ^ ' 



where the matrix Xi depends on the occupation r.; of site i 

X, = nD + (1 - n)E (49) 

and the matrices D and E satisfy the following algebraic rules 

DE- ED = D + E 

{W\{aE --/D) = {W\ (50) 
{f3D - SE)\V) = \V) . 
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Let us check on the simple example of figure[8]that expression (|48l) does give the steady state weights: if one chooses 
the configuration where the first p sites on the left are occupied and the remaining L — p sites on the right are empty, 
the weight of this configuration is given by 

{W\DPE^-P\V) 

{W\(D + E)i^\V) ■ ^ ' 

For (|48)) to be the weights of aU configurations in the steady steady, one needs that the rate at which the system 



I I I \ \ I I 

Figure 8: The three configurations which appear on the left hand side of (|52p and from which one can jump to the configuration 
which appears on the right hand side of (I52|l . 

enters each configuration and the rate at which the system leaves it should be equal. In the case of the configuration 
whose weight is (|51|) . this means that the following steady state identity should be satisfied (see figure [8]): 



{W\EDP-^E^-P\V) {W\DP-^EDE^-P-^\V) {W\DPE^-p-'^D\V) _ {W\DPE^-p\V) 
" {W\iD + E)^\V) ^ {W\{D + E)^\V) ^ ^^WliDTWWr ~ ' {W\{D + E)i^\V) 

This equality is easy to check by rewriting ([5^ as 

{W\{aE - jD)DP-^E^~P\V) {W\DP-^{DE - ED)E^-p-^\V) (W\DPE^-p-\i3D - 5E)\V) 

{wKdTWW) (wKdTWW) {W\{D + E)^\V) ~ 



(53) 



and by using ([50)1 . A similar reasoning [10| allows one to prove that the corresponding steady state identity holds for 
any other configuration. 

A priori one should construct the matrices D and E (which might be infinite-dimensional [13]) and the vectors {W\ 
and \V) satisfying ([50]) to calculate the weights ((48|) of the microscopic configurations. However these weights do not 
depend on the particular representation chosen and can be calculated directly from (|50p. This can be easily seen by 
using the two matrices A and B defined by 

A = PD-SE 

B = aE-jD (54) 

which satisfy 

AB -BA^ {a(3 - jS){D + E) = {a + S)A + (/? + -f)B . (55) 



Each product of D's and i?'s can be written as a sum of products of A's and B's which can be ordered using ([SS]) 
by pushing all the A's to the right and all the B's to the left. One then gets a sum of terms of the form BPA'^, the 
matrix elements of which can be evaluated easily ((VF|i3PA'|y) = (VF|y)) from ([50)1 and ([5l 
One can calculate with the weights the average density profile 

_ {W\{D + Ey-^D{D + E)^-'\V) 
~ {W\{D + E)L\V) 

as well as all the correlation functions 

{W\{D + Ef-'^D(D + Ey-'-^D{D + E)^-^V) 



{W\iD + E)L\V) 
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and one can recover that way pTj) and ([44]). 

Using the fact that the average current between sites i and i + 1 is given by 

_ {W\{D + Ey-\DE - ED){D + E)^-'-^\V) _ {W\{D + E)^-^\V) 
^'^ {W\iD + E)L\V) ~ {W\{D + E)i^\V) 

(of course in the steady state the current does not depend on i) and from the expression (j43p one can calculate the 
normalization 

{W\iD + E)^\V) ^ 1 + 

{w\v) iPa-Pt)^ r(^ + ^) ^ > 

(see equation (3.11) of 26] for an alternative derivation of this expression). 

Remark: when Pa = Pt = the two reservoirs are at the same density and the steady state becomes the 
equilibrium (Gibbs) state of the lattice gas at this density r. In this case, the weights of the configurations are those 
of a Bernoulli measure at density r, that is 

L 

Pro({Ti,T2,...ri}) = [][rr, + (l-r)(l-r,)] . (57) 

i=l 

This case corresponds to a limit where the matrices D and E commute (it can be recovered by making all the 
calculations with the matrices (|48l50p for pa ^ Ph and by taking the limit pa Pb in the final expressions, as all the 
expectations, for a lattice of finite size L, are rational fimctions of pa and pb). 



VII. ADDITIVITY AS A CONSEQUENCE OF THE MATRIX ANSATZ 

As in (|48[) the weight of each configuration is written as the matrix element of a product of L matrices, one can 
try to insert at a position Li a complete basis in order to relate the properties of a lattice of L sites to those of two 
subsystems of sizes Li and L — Li. To do so let us define the following left and right eigenvectors of the operators 
PaE - (1 - pa)D and (1 - pb)D - p^E 

{pa,a\ [paE - (1 - Pa)D] = a{pa,a\ 

[(1 - pb)D - pbE] \pb, b) = b\pb, b) . (58) 
It is easy to see, using the definition 0, that the vectors {W\ and \V) are given by 

{W\ = (p,,(a + 7)-i| 

\V)^\pb,{P + 5)-^) ■ (59) 

It is then possible to show, using simply the fact l[50|) that DE — ED — D + E and the definition of the eigenvectors 
([55)) . that (for pb < Pa) 

{pa,a\Y,Y2\pb,b) _ I dp {Pa-PbY+'' {pa,a\Y^p,b) {p,l-b\Y2\pb,b) ^gp^ 



{pa,a\pb,b) J 2iTT {pa- pY+''{p- Pb) {pa,a\p,b) {p,l-b\pb,b) 

Pb<\p\<Pa 

where Yi and Y2 are arbitrary polynomials of the matrices D and E. 

Proof of dM]): to prove it is suflicient to choose Yi of the form [paE - {I - pa)D]"[D + E]''' and Y2 

of the form [D + i?]" [(1 — pb)D — Ph-E]" (one can show, using DE — ED = D + E, that any polynomial Yi or Y2 
can be reduced to a finite sum of such terms). Then proving ([60)l for such choices of Yi and Y2 reduces to proving 

{pa,a\iD + E)"'+""\pb,b) / dp {pa~Pb)''+'' {pa,a\{D + E)"'\p,b) {p,l-b\{D + E)""\pb,b) 

- ^ - . (61) 



{pa,a\pb,b) J 2iTr {pa~ p)''+''{p- Pb) {pa,a\p,b) {p,l-b\pb,b) 

Pb<\p\<Pa 
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As from (ISH)) one has 

{pa,a\{D + E)^\p,b) _ T{L + a + b) 



{pa,a\p,b) {pa- Pb)^ ^{a + b) ' 

Then (f60|) and (I61|) follow as one can easily check that 

r{n' + n" + a + b) f dp {pa - PbT^'' r{n' + a + b)T{n" + 1) 



{pa - Pby''+"" J 2iTT {pa - {p - pi,}"" + ^ T{a + b) 

Pb<\p\<Pa 

If one normalizes ([60|) by (l56|) one gets 

(Pa,a|yir2|pb,6> r(L + L' + a + 6) / dp (pa - 



(62) 



(63) 



{pa,a\{D + E)L+^'\pb,b) r{L + a + b)r{L' + 1) J 2nT [p a - pY+''+'^{p - pbY+^' 

Pb<\p\<Pa 

(pa,a|Yi|p,5) {p,l-b\Y2\pb,b) 



(64) 



(Pa,a|(i? + i?)^|p,fe> (p,l-6|(i? + i?)^'|pb,6) ■ 

An additivity relation more general than (pO]) can be proved for the ASEP 53]. The special case of the TASEP 
will be discussed in section IXVII below. 

VIII. LARGE DEVIATION FUNCTION OF DENSITY PROFILES 

If one divides a chain of L sites into n boxes of linear size I (there are of course n = L/l such boxes), one can try to 
determine the probability of finding a certain density profile {pi, p2, ■■■Pn}, i-e. the probability of seeing Ipi particles 
in the first box, lp2 particles in the second box, ... ^p„ in the nth box. For large L one expects (fT5|) the following L 
dependence of this probability 

PrOL{pi,...Pn\pa,Pb) ^ exphLTn(pi,P2,--Pn\Pa,Pb)] ■ (65) 

If one defines a reduced coordinate x by 

i — Lx (66) 

and if one takes the limit I — > oo with I <^ L so that the number of boxes becomes infinite, one gets as in ([T7)) the 
large deviation functional J-{p(x)) 

PrOi({p(x)}) ^ exp[-LT{{p{x)}\pa,Pb)] ■ (67) 
For the SSEP (in one dimension), the functional J-{p{x)\pa, Pb) is given by the following exact expressions: 
at equilibrium, i.e. ior pa = Pb = r 

T{{p{x)}\r,r)= l\{p{x),r)dx (68) 

where 







B{p,r) = {l-p) log ^ + p log ^ . (69) 
1 — r r 

This can be derived easily. When Pa = Pb = r, the steady state is a Bernoulli measure ([57]) where all the sites are 
occupied independently with probability r. Therefore if one divides a chain of length L into L/l intervals of length 
one has 
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and using Stirling's formula one gets (|68l69p . 

For the non-equilibrium case, i.e. for pa 7^ Pb, it was shown in [2^, [2^, [S^l that 

F'{x) 



^{{p{x)}\Pa, Pb) = / dx 

Jo 



B{pix),Fix)) + \og 



Pb - Pa 



where the function F{x) is the monotone solution of the differential equation 

F(l - F)F" 



p{x) = F + 



pi2 



(71) 



(72) 



satisfying the boundary conditions F(0) = pa and F(l) = pf,. This expression shows that is a non-local functional 
of the density profile p{x) as F(x) depends on the profile p{y) at all points y. For example if the difference pa — Pb 
is small, one can expand ^ and obtain the expression ([27|) where the non-local character of the functional is clearly 
visible: at second order in Pa — Pb, one gets by solving ([72|) 



F ^Pa - [Pa - Pb)x 



{Pa ~ Pbf 
Pa(l - Pa) 



--p*{x) - 



{Pa - Pb)^ 
Pa(l - Pa) 



(l-x) I y{p{y) ~ pa)dy + X I {1 - y){p{y) - pa)dy 

1 



(l-x) y{p{y) - p*{x))dy + X {I - y){p{y) ~ p* {x))dy 



0((Pa-P6)') 

0{{pa~Pb)^) (73) 



and this leads to ((27|l by replacing into ((7T|) . 



Derivation of (j71ll72p : in the original derivation of (|71l72p from the matrix ansatz |25l |26 l| the idea was to 
decompose the chain into L jl boxes of I sites and to sum the weights given by the matrix ansatz (|48l50p over all the 
microscopic configurations for which the number of particles is lp\ in the first box, lp2 in the second box Ip^ in 
the nth box. 



An easier way of deriving (|71l72p is to write (we do it here in the particular case where a + 6 = 1, i.e. ^ 



+7 ' /3+<5 



1, 



and Pb < Pa) from (j60p and (j56p where Yi and I2 represent sums over all configurations with kl sites with a density 
pi in the first I sites, ... pk in the fc-th I sites, and Y2 a similar sum for the (n — k)l remaining sites. 



Pnlipi 

, P2---Pn \Pa 1 Pb) 



(kiy. {{n~k)l)\ 



{nl)\ 



dp 
2i7r' 



(74) 



Pb<\p\<Pa 



l^p p^yil+1 

(Pl • ■ -Pfc |Pa , P) Pin-k)l (Pk+l .■■Pn\p,Pb) 

Note that in ([71]) the density p has become a complex variable. This is not a difficulty as all the weighs (and therefore 
the probabilities which appear in (j74p ) are rational functions of pa and pb. 

For large n/, if one writes k = nx, one gets by evaluating (|74p at the saddle point 



^n{Pl:P2, ■■■Pn\Pa,Pb) = Hiax xJ^kiPl, ■■■Pk\Pa,F) + (1 - x)J^n-kiPk+l, ■■■Pn\F,Pb) 

Pb<F<Pa 

+ x\og (^^) + (1 - ^) log i^Y^) - log(P« - Pb) 



(75) 



(To estimate (1741) by a saddle point method, one should find the value F oi p which maximizes the integrand over 
the contour. As the contour is perpendicular to the real axis at their crossing point, this becomes a minimum when 
p varies along the real axis). If one repeats the same procedure n times, one gets 



J^n{pi,P2T--Pn\Pa,Pb) ^ Hiax 

Pi = Fo<Fi..<Fi<..<F^=p^ n 



^^Ti IP {Fr^i-Fi)n \ 

-y^^i {Pi\Fi-i, F,) + log 

nf^l \ pa-pb J 



(76) 



For large n, as Fi is monotone, the difference — Fi is small for almost all i and one can replace J-'i{pi\Fi-i, Fi) 
by its equilibriimr value J-'i{pi\Fi, Fi) = B {pi , Fi) . If one write Fi as a function of i/n 



F,^F\- 

n 



(77) 
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d76l) becomes 



J^{{p{x)}\pa, Pb) = max / dx 

F{^) Jo 



B{p{x),F{x))+\og 



F'{x) 



(78) 



Pb- Pa. 

where the maximun is over all the monotone fmictions F{x) which satisfy F{Q — pa and F{\) — pb and one gets (|71l72p . 



Remark: One can easily get from (|71|72p the generating function Q{{a{x)}) of the density ((30|) for the 
SSEP: 



Q{{a{x)}) = / dx 
Iq 



log(l-F + i^e"(^)) - log- 



F' 



Pb - Pa 



where F is the monotone solution of 



F 







1 - F + Fe"(^) 

with F{0) = Pa and F{1) = pb- For small a{x) the solution of ([80)1 is to second order in the difference pa — Pb 



= P*{x) - {pa - Pbf 



(l-x) I ya{y)dy + x\ (l-y)a{y)dy 



(79) 



(80) 



(81) 



This leads to Q{a{x)) at order {pa — pb) 



g{{a{x)}) = / dx 



p*{x)a{x) + 



p*{x){l~p*{x)) 



alxY 



- {pa - Pb) I dx I dy x{l - y) a{x)a{y) (82) 



and one recovers through ([33]) the expression of the two-point correlation function ([29 



IX. THE MACROSCOPIC FLUCTUATION THEORY 



For a general diffusive one dimensional system (figure [2]) of linear size L the average current and the fluctuations 
of this current near equilibrium can be characterized by two quantities D{p) and cr(p) defined by 



lim 

t— >CXD t 



{Qt) _ D{p) 



L 



[Pa - Pb) for {pa - Pb) small 



(83) 



limM) 

t— >oo t 



<y{p) 



for Pa = Pb 



(84) 



where Qt is the total number of particles transferred from the left reservoir to the system during time t. 

Starting from the hydrodynamic large deviation theory [2ll, [t^I Bertini, De Sole, Gabrielli, Jona-Lasinio and 
Landim [4ll . |4^ . [4^ have developed a general approach, the macroscopic fluctuation theory, to calculate the large 
deviation functional T of the density (jl7p in the non-equilibrium steady state of a system in contact with two (or 
more) reservoirs as in figure [21 Let us briefly sketch their approach. For diffusive systems (such as the SSEP), the 
density Piit) near position i at time t and the total flux Qi{t) flowing through position i between time and time t 
are for a large system of size L and for times of order L^, scaling functions of the form 



P.{t) 



i t 



and 



,{t) = LQ 



i t 



(85) 



(Note that due to the conservation of the number of particles, the scaling form of pi{t) implies the scaling form of 
Qi{t)). If one introduces the instantaneous (rescaled) current defined by 



3{x, t) 



dQ{x, t) 



(86) 
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the conservation of the number of particles imphes that 



dp{x,T) d'^Q{x,T) _ dj{x,T) 



dr 



drdx 



dx 



(87) 



Note that the total flux of particles through position i = [La;] during the macroscopic time interval dr, i.e. during the 
microscopic time interval L^dr, is Lj(x,T)dT. Thus the microscopic current is of order 1/L while the rescaled current 
j remains of order 1. 

The macroscopic fluctuation theory [4ll [4^ l43j starts from the probability of observing a certain density profile 
p{x, t) and current profile j (x, r) over the rescaled time interval ti < t < T2 



Qri.T2(^{pix,T),j{x,T)} 



exp 



-L [ dr' [ dx 
Jti Jo 



j(x,r') + D{p{x,r'))^ 



/\\d£(xy) 



2a{p{x,T')) 



where the current j{x, s) is related to the density profile p{x, s) by the conservation law (I87p and the functions D{p) 
and cr(p) are defined by ()83|84p . Note that a similar expression was obtained in |82] by considering stochastic 
models in the context of shot noise in mesoscopic quantum conductors. 

Then Bertini et al [41] show that to calculate the probability of observing a density profile p{x) in the steady state, 
one has to find out how this deviation is produced. For large L, one has to find the optimal path {p{x, s),j{x, s)} for 
— oo < s < T in the space of density and current profiles and 

(89) 



(90) 



Pro(p(x)) max Q_oo,r ( {p(a;, s), i(a;, s)} 

{p(2:,s) J(2:,s)} ^ 

which goes from the typical profile p* (x) to the desired profile 

p{x,-oo) = p*{x) ; p{x,t) = p{x) . 
This means that the large deviation functional J- of the density ([67]) is given by 



J-{p(x)) = min / dr' I dx- 

{p(x,s),j(x.s)} J -oo Jo 



j{x,T')+D{pix,T')) 



dx 



2a{p{x,T>)) 



(91) 



where the density and the current profiles satisfy the conservation law (|87p and the boundary conditions ([50)1 . 

Finding this optimal path p(x,s),j{x,s) with the boundary conditions (|90p is usually a hard problem. Bertini et 
al [4l| were however able to write an equation satisfied by J^: as ((9T|) does not depend on r, one can rewrite it as 



T{p{x)) = min 



2<T(f(l)) 



where p{x) — Sp{x) — p{x, t — dr) and j{x) = j{x, r). Then if one defines U{x) by 

dT{{p{x)}) 



U{x) 



5p{x) 



(92) 



(93) 



and one uses the conservation law Sp{x) 
is given by 



dx 



dr one should have according to ([92]) that the optimal current j{x) 



jix) = -D{pix))p'ix) + <j{p{x))U'{x) 
Therefore starting with p(x, t) — p{x) and using the time evolution 

d'p{x,s) dj{x,s) 



ds 



dx 



(94) 



(95) 



with j related to p by ([M]) one should get the whole time dependent optimal profile p{x, s) which converges to p*{x) 
in the limit s — s- — oo. The problem of course is that U{x) defined in (j93p is not known. One way of thinking of 
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the problem is to say that for every p{x), the function U{x) should be adjusted so that the dynamics (|95l94p gives 
p{x, s) — > p*{x) as s ^ — oo. 

One can write from (|92p (after an integration by part and using the fact that C/(0) = C/(l) = if p(0) — pa and 
p{l) = ph) the equation satisfied by U'{x) 



f 

Jo 



dx 



Dp' 



-U'\ 



Dp' 



1-^ 



(96) 



which is the Hamilton- Jacobi [41j equation of Bertini et al. For general D{p) and (j{p) one does not know how to find 
the solution U'{x) of (|96p for an arbitrary p{x) and thus one does not know how to get a more explicit expression of 
the large deviation function T{{p{x)}). 

One can however check rather easily whether a given expression of !F{{p(x)}) satisfies (|M)) since U'(x) can be 
calculated from For the SSEP one gets from (|78l93p 



U{x) = log 



p{x){l~F{x)) 



{l^p{x))Fix) 

is indeed satisfied using the expressions of D 



(97) 

= 1 



with F{x) related to p{x) by ([7^ . One can then check that 
and cr = 2p(l - p) for the SSEP (see (fTTC)) below). 

In fact J- is known, one can obtain the whole optimal path p{x, s) from the evolution (j95p with j related to p by 
([Ml which becomes for the SSEP 



dp(x, s) 
dx 



+ (T(p(a;,s))log 



p{x,s){l~F{x,s)) 
{\-p{x,s))F{x,s) 



(98) 



where F is related to p by ([7^ . For (|71l72p to coincide with ([?!]) , the optimal profile p evolving according to 
should converge to p*{x) as s — ^ — oo. One can check that this evolution of p[x,s) is equivalent to the following 
evolution ([H) of F 



dF{x,s) _ d^F{x,s) 



ds 



ds 



(99) 



where F is related to p by ([7T|) . Clearly F{x,s) — + p*{x) and therefore p{x,s) P*{x) as s ^ — oo. Thus (|95l98p 
do give the optimal path in (PT|) with the right boundary conditions ([M)) and ([OT]) coincides for the SSEP with the 
prediction (|71l72p of the matrix approach. 

Apart for the SSEP, the larg e deviation functional J-' of the density is so far known only in very few cases: the Kipnis 
Marchioro Presutti model j83l. [83| the weakly asymmetric exlcusion process [13, and the ABC model [bI, [T^I on a 
ring. 



X. LARGE DEVIATION OF THE CURRENT AND THE FLUCTUATION THEOREM 

For a system in contact with two reservoirs at densities pa and ph, as in figure[21 one can try to study the probability 
distribution (86j of the total number Qt of particles which flows through the system during time t. For finite t, this 
distribution depends on the initial condition of the system as well as on the place where the flux Qt is measured 
(along an arbitrary section of the system, at the boundary with the left reservoir or at the boundary with the right 
reservoir). In the long time limit however, if the system has a finite relaxation time and if the number of particles in 
the system is bounded (i.e. infinitely many particles cannot accumulate in the system) the probability distribution of 
Qt takes the form 




(100) 



where the large deviation function F{j) of the current j depends neither on the initial condition nor on where the flux 
Qt is measured. This large deviation function F{j) has usually a shape similar to ap{p) in flgurelHl with a minimum 
the typical value j* = (J) (the avergae current) where F{j*) = 0. 

It is often as convenient to work with the generating function (e^^^Y the long time limit 

(e^Q') ~ e^(^) * (101) 
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where /i(A) is clearly the Legendre transform of the large deviation function F(j) 

fiiX) = max[Aj - F{j)] 

As in section HVl . the knowledge of //(A) determines the cumulants of Qt 



lim 

t—*oo 



(102) 



(103) 



A=0 



when the expansion in powers of A is justified. 



l89l . l9d . l9l| the large deviation func- 



According to the fluctuation theorem, [45|, |46|, |47|, |48|, |49|, [STI, |58|, \m, 
tion F{j) of the current satisfies the following symmetry property 

FU) - Fi-j) ^ -j[^og Za- log Zb] and ^(A) = /i (-A + logz,, - logz^) . (104) 

Proof: Following previous derivations [i^ [s^ for stochastic dynamics the fluctuation theorem (|104p can be easily 
recovered [sgI ) from the generalized detailed balance relation ([9]). This can be seen by comparing the probabilities of 
a trajectory in phase space and of its time reversal for a system in contact with two reservoirs. A trajectory ^^Traf^ 
is specified by a sequence of successive configurations Ci, ...Ck visited by the system, the times ii, ...tk spent in each 
of these configurations, and the number of particles qa.i,Qb,i transferred from the reservoirs to the system when the 
system jumps from Ci to C^+i. 



Pro(Traj) = dt 



k~l 



k-l 



exp 



Y^t, ria 



where r(C) = J2c" q,, ^qa,qbi^' ^ ^) ^'^d dt is the infinitesimal time interval over which jumps occur. 

For the trajectory " —Traj" obtained from ^^Traj" by time reversal, i.e. for which the system visits successively the 
configurations Ck, ■■■Ci, exchanging —qa,i,—Clb,i particles with the reservoirs each time the system jumps from C^+i 
to Ci, one has 



Pro(-rmj) = dt 



k-l 



■fe-1 



Y[w. 



-96. i (Cjj C'i+l) 



exp 



One can see from the generalized detailed balance relation ([9]) that 

exp 



Pro(rraj') 



where Qj 

the system during time t 



PTo{—Traj) 

(b) _ 



QI"^ log Za - Q[''^ log Zb 



(105) 



qa,i log Za - Qb.i log Zb = CXp 

J2i Qa,i and Q\"' — J2i 1b,i are the total number of particles transferred from the reservoirs a and b to 
ing time t. 

In general Qt""* and Qj^'' grow with time but their sum remains bounded (if one assumes that particles cannot 



accumulate in the system). Therefore for large time Qt = Q'f^ — — qI''' 



PTo{Traj) 



o{t) and 
exp [Qt{logZa - log Zb)] 



(106) 



Pro(— Traj) 

Summing over all trajectories [5^, taking the log and then the long time limit (jlOOp leads to the fluctuation theorem 
(fT04l) . 

Remark: The fluctuation theorem predicts a symmetry relation similar to (jl04p for the heat current for a 
system in contact with two heat baths at unequal temperatures as in figure [TJ Under similar conditions as for the 
current of particles (the energy of the system is bounded and the relaxation time is finite - see [q^, [l^ 111] for 
counter-examples where the energy is not bounded in which case the fluctuation theorem has to be modified) one 
gets that the distribution of the energy Qt flowing through the system during a long time t is given by (jlOOp and 
that the large deviation function F(j) or its Legendre transform satisfy 



Fij) - Fi-j) = -J 



1 



1 



/i(A) =M -A + 



1 



(107) 



which states that the difference F{j) — F{—j) is linear in j with a universal slope related to the difference of the 
inverse temperatures. Note that — pf-) is the rate of entropy production which is the quantity generally used 



to state the fluctuation theorem [43, WQl ^ 
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XI. THE FLUCTUATION-DISSIPATION THEOREM 



In the limit of small Tq — Tf, (i.e. close to equilibrium), one can recover from p07|) the fluctuation- dissipation 
relation between the response to a small temperature gradient 

^ ^ (T, ~ n)D for T, - n small (108) 
and the variance of the energy flux at equilibrium 

Ml 
t 

In fact from these definitions of D and a, one has 



a for Ta=Tb. (109) 



fiiX) = (Ta - n)DX + + O (A3, A2(T, - T,), A(T, - Tf,)^) (110) 

and using the fluctuation theorem (|107p . one gets that the coefficients a and D have to satisfy 

a = 2kT^D (111) 
which is the usual fluctuation-dissipation relation. In general both D and a depend on the temperature T^. 
The same close-to-equilibrium expansion of (|104p for a current of particles leads to 

a = = 2kTp^K(p)D (112) 

alogz 

where the coefficients a and D are defined as in (|109ll08p by 



(Q*> ... „^Fs „ „_.,^„u (Q' 



{Pa - Pb)D for Pa - Pb smaU and — > cr for pa = (113) 

I- t 

and where k(/9) is the compressibility (|26|) at equilibrium. 

To see why the compressibility appears in (|112p . one can write logZ = —F/kT = —Vf{N/V)/kT 
where F is the total free energy and / the free energy per unit volume. One then uses the facts that 
the fugacity z is given by fcTlogz — dF/dN — f'{p), that p = —dF/dV = pf'{p) — f{p) and thus 
dp/d\ogz = kT/f"{p) = kTpdp/dp ^ kTp^K{p). 

In the case of the SSEP, one has from (|43l83p that Z^ssep = 1- As the free energy f{p) of the SSEP (at 
equilibrium at density p) is 

f{p) = kT[p\ogp+ (1 - p) log(l - p)] (114) 

one has 

logz = log—^ . (115) 
1 - p 

Thus dp/d\ogz = kT/f"{p) = p{l - p) and thanks to fm]) and (|43l83p one gets 

i'ssEP = 1 ; crssEP = 2/9(1 - p) (116) 



Note that in (|83l84p there is, compared with (|109ll08lll3p . an extra \/L factor in the definition of a and D to get a 
finite large L limit of <j and D. Of course, with this extra 1/L fator, both piip and (|112p remain valid. 
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XII. CURRENT FLUCTUATIONS IN THE SSEP 



For the SSEP, if Qt is the total number of particles transferred from the left reservoir to the system during a long 
time i, one has ()101|) 

(e^Q*) ^e'^^^)* . (117) 
The fluctuation theorem (jl04p implies ()115|) a symmetry relation satisfied by ^(A) 

MA)=A^f-A-log-^+log-^) . (118) 

but of course this symmetry does not determine p(A). We are now going to see that, because the evolution is 
Markovian, /i(A) can be determined as the largest eigenvalue of a certain matrix [sol [osl [qgHqtI. 
The probability Pt{C) of finding the system in a configuration C at time t evolves according to (jU 



dPtjC) 
dt 



J2 W{C, C')Pt{C') - W{C', C)Pt{C) . (119) 



Among all the matrix elements W{C, C"), some correspond to exchanges of particles with the left reservoir and others 
represent internal moves in the bulk or exchanges with the right reservoir. Thus one can decompose the matrix 
W{C, C) into three matrices 

WiC,C')^WiiC,C') + WoiC,C') + W^iiC,C') (120) 

where here the index is the number of particles transferred from the left reservoir to the system during time dt, when 
the system jumps from the configuration C to the configuration C. One can then show [50, 95, 96] that /i(A) is 
simply the largest eigenvalue (more precisely the eigenvalue with largest real part) of the matrix Mx defined by 

Ma(C, C") = e^WiiC, C) + Wo{C, C) + e-^W-i{C, C) - S{C, C) ^ W{C", C) . (121) 

C" 

In fact the joint probability Pt{C, Qt) of C and Qt evolves acording to 

dPtiC^Qt) ^ W,iC,C')PtiC\Qt-q)-J2wiC\C)PtiC,Qt). (122) 

C ij=-l,0,l c 

Then if Pt{C) = e-^<3'Pf (C, Qt) one has 

^Iji^lRA ^J2Mx{C,C")PtiC') . (123) 
and this shows that /i(A) is the eigenvalue with largest real part of the matrix M\. 

The size of the matrix M\ grows like 2^ (which is the total number of possible configurations of a chain of 
L sites). In |53| a pertubative approach was developed to calculate /i(A) in powers of A. Let us sketch briefly this 
approach: one can write down exact expressions for the time evolution ^e'^*^*) or of (^e^^* H{C)') where H{C) is an 
arbitrary function of the conflguration C at time t. For example 

Qt with probability 1 ^ a{\ — Ti)dt — ^Tidt 
Qt+dt = "( Qt + 1 with probability a(l - Ti)dt (124) 
Qt — 1 with probability ^Tidt 



and therefore 



^^i^=a(e^-l)((l-ri)e^«')+7(e-'-l)(ne^'5') • (125) 
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Similarly one can show that for 1 < i < L 

^^^^^ = a{e' ^ 1) ((1 - r,)ne'Q^)+j{e-' ~ 1) (nne'Q^) + ((r.+i - 2t, + r,_i)e^Q') (126) 
the cases i = 1 are i = L being slightly different 



dt 



ae 



((1 - Ti)e^Q') - 7(Tie^Q') + ((r2 - ri)e^Q') (127) 



^^^^^ = a(e^-l) ((1 - r^)rLe'Q')-j{e'^-l) {nri^e^Q^) + {{ri^^, - ri)e^Q') +5 ((1 - ri)e^Q*)-/3(r,.e^«0 . 

(128) 

In the long time limit (e'^'^*) ~ e''^'*'''* and one can define a measure {.)\ on the configurations C 



{H{C))x = hm ^ )' ' (129) 



From p25m28p one gets 

/i(A) = «(e^ - 1) ((1 - n))^ + 7(e-^ - 1) (n), (130) 

^(A)(t,)a = a(e^ - 1) ((1 - ri)T,);, + 7(6"^ - 1) (nr,);, + ((t,+i - 2r, + t,_i))^ (131) 

Ai(A)(ri)A = ae^(l - ti))a - 7 (n), + ((r^ - n)), (132) 



MA)(rL)A - «(e' - 1) ((1 - T^)TL)y^ - 7(e-^ - I) (rir^)^ + ((r^^i - r^))^ + 5 ((1 - r^));, - /3 (r^)^ 



(133) 



We see that to get /i(A) at order A*^, one needs to know (I130p the one-point function (Ti)\ at order A*"'^^, the two 
point fmictions {TiTj)\ at order A*^"^ (see (|13mi33|) ) and so on up to the fc— point functions at order A°. As the steady 
state weights P{C) for the SSEP are known exactly (section 25]) [l^ilH, the correlation functions are known 

at order A" and one can truncate the hierarchy at the level of the fc— point functions. 

In [s^l this perturbation theory based on the hierarchy (|130m33|) was developed to calculate ^{\) in powers of A. 
The main outcome of this perturbation theory |50| is that /i(A), which in principle depends on L, A and on the four 
parameters a, (3, 7, (5, takes for large L a simple form 



A*(A)--i?(^)+0( — ) (134) 



- (e^ - l)pa + [e-^ - l)pb - (e^ - 1)(1 - e-^)paPb ■ (135) 



where uj is defined by 

where pa and ph are given in ([2|). The perturbation theory gives up to fourth order in uj 

i?H=--y + ^-|^ + 0(a.^). (136) 

The fact that p{\) depends only on p^j Pb and A through the single parameter lo is the outcome of the calculation, 
but so far there is no physical explanation why it is so. However uo remains unchanged under a number of symmetries 
[50| (left-right, particle-hole, the Gallavotti-Cohen symmetry (IllSp *) implying that p{\) remains unchanged as it 
should under these symmetries. 

From the knowledge of R{oj) up to fourth order in lo, one can determine ^] the first four cumulants (|103p of the 
integrated current Qt for arbitrary pa and pb'. 
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• For Pa = 1 and pf, = 0, one finds 



{Qt 



{QT, 



m 
t 



_ 1 

~ T 
1 

" 3l 
1 

Tsl 

-1 



o 



o 



o 



105i 



o 



(137) 
(138) 
(139) 
(140) 



These cumulants are the same as the ones known for a different problem of current flow: the case of non- 
interacting fermions through a mesoscopic disordered conductor |98l. 1991 . This can be understood as a theory 
[S^ similar to the macroscopic fluctuation theory of Bertini et al [4l|, [i^, IH, can be written for these 
mesoscopic conductors with the same D{p) and a{p) as for the SSEP l|116p . 

• For Pa — Pb = \ which corresponds to an equilibrium case with the same density 1/2 in the two reservoirs, one 
finds that all odd cumulants vanish as they should and that 



(Q 



1 

2L 



O 



{Q 



t 



o 



1 

1 

L2 



(141) 
(142) 



Because p{X) depends on the parameters pa, Pb and A through the single parameter uj, if one knows p{X) for one single 
choice of pa and pt, then p35ll36p determine p{X) for all other choices of pa,Pb- In 50], it was conjectured that for 
the particular case Pa ^ Pb — ^, not only the fourth cumulant vanishes as in p42p . but also all the higher cumulants 
vanish, so that the distribution of Qt is Gaussian (to leading order in 1/L). This fully determines the function R{llj) 
to be 



R{lu) = [\og{VT 



(143) 



One can then check that, with this expression of R{u!), not only p37ll40p but all the higher cumulants of Qt in the 
case Pa = 1 and ph ^ coincide with those of fermions through mesoscopic conductors [sol [98| . 



XIII. THE ADDITIVITY PRINCIPLE 

In (5]| . another conjecture, the additivity principle based on a simpler physical interpretation, was formulated 
which leads for the SSEP to the same expression (|135ll43p as predicted in section IXIII and can be generalized to 
obtain F[j) or /Lt(A) for more general diffusive systems. 




Figure 9: In the additivity principle one tries to relate the large deviation function of the current Fl+l' (j) of a system of size 
L + L' to the large deviation functions F^ij) and Fi^i{j) of its two subsystems (|146|) . 
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For a system of length L + L' in contact with two reservoirs of particles at densities pa and pb, the probability of 
observing, during a long time t, an integrated current Qt = jt has the following form (jlOOp 

PrOL+L' (j, Pa, Pb) - e-*^-+-'(j"'''-^'') . (144) 

The idea of the additivity principle is to try to relate the large deviation function F^+l' (j, Pa, Pb) of the current to 
the large deviation functions of subsystems by writing that for large t 

ProL+L' U, Pa, Pb) ^ max [ProL (j, Pa, p) x PrOi' (j, p, pb)] . (145) 
p 

This means that the probability of transporting a current j over a distance L + L' between two reservoirs at densities 
Pa and Pb is the same (up to boundary effects which give for large L subleading contributions) as the probability of 
transporting the same current j over a distance L between two reservoirs at densities pa and p times the probability 
of transporting the current j over a distance L' between two reservoirs at densities p and pb. One can then argue 
that choosing the optimal p makes this probability maximum. From (|145p one gets the following additivity property 
of the large deviation function 

Fl+L' {j, Pa, Pb) = max [Fl (j, Pa, p) + Fl' {j, p, Pb)] ■ (146) 
p 

Suppose that we consider a diffusive system for which we know the two functions D{p) and a{p) defined in (j83l84p . 
If one accepts the additivity property (jl46[) of the large deviation function, one can cut the system into more and 
more pieces so that 



'n-l 



Fl (j, Pa, Pb) ^ min < V FL/„{j, Pt,Pt+i) 



(147) 



. i=0 



where po = pa and p„ = p&. 

For large n, the optimal choice of the p^'s is such that the differences pi — pi+i are small. For a current j of order 
1/L, and for pi — p^+i small, one can replace ()83|84p Fi for I = L/k hy 



Fi{3,Pi,Pt+i) 



-P(P. )(Pi-Pi + l) j2 



lipA 



(148) 



and by taking the limit n ^ oo (keeping I = L/n large for (|83l84p to be still valid) one gets [5l| 



Fl {j,Pa,Pb) = y max 

1j p{x) 



' [jL + D{p)p'- 

Mp) 



-dx 



(149) 



where the optimal profile p{x) (for large n, the optimal pi in p47p ) given by pi — p(i/n) should satisfy p(0) = pa and 
P(l) = Pfc- 

One can show |5ll . |56| that the optimal profile in (jl49p (when pa ^ Pb and the deviation of current j is small enough 
for this optimal profile to be still monotone) is given by 



iLj)^l + 2Ka{p{x)) 



p'{x)^ = 



(150) 



D^{p{x)) 

where the constant K is adjusted to insure that p(0) = pa and p{l) = pb- Replacing p{x) by (|150p in (|149p leads to 



FL{j,Pa,Pb) = j 



1 + Ka{p) 
[l + 2i^a(p)]i/2 



1 



Dip) 
<^{p) 



dp 



where the constant K is fixed from (|150p by the boundary conditions (p(0) — pa and p(l) — pb) 

r"' Dip) 



Pb 



[l + 2Xa(p)]i/ 



(151) 



(152) 



Expressions p51ll52p give therefore F^ij) in a parametric form. 
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The optimal profile (jl49p remains unchanged when j 
(|151|152p and one gets that (fm]) 



-j (simply the sign of [1 + 2fCcr(/3)]^/^ is changed) in 



Pb 



dp = -j(log Za - log Zb) 



(153) 



Thus the expression (|151ll52p does satisfy the fluctuation theorem (|104p . 

From (|151ll52p one can calculate /Lt(A) by (|102p and one gets [5l|, |5f| a parametric form 



M(A,Pa,Pb 

with K = K{X, Pa, Pb) is the solution of 

A = 



K 



D{P) dp 
VI + 2Kcj{p) 



(154) 



Pb 



Vl + 2Ka{p) 



(155) 



One can then get [51| . by eliminating K (perturbatively in A) the expansion of /i(A) in powers of A and therefore the 
cumulants (|103p in the long time limit for arbitrary pa and pb 



{Qt 



where 



L^^' 1 



ill 



{Q 



^3(/3/i-/|) {Qt 
L 



If 



D[p) a{p)"-' dp 



t 



1 3(5/4/1 
L 



14/1/2/3 



9/1) 



I! 



(156) 



(157) 



Using the fact ([TTC)) that for the SSEP D{p) = 1 and a{p) = 2p(l - p), one can recover [5l| from (|154ll55ll56p the 
abov e expressions (|137m40ll43l) . The validity of (|151ll52p has also been checked for weakly interacting lattice gases 

[Tool. 



We have seen in IIXI that the macroscopic fluctuation theory gives the probability ([55)1 of arbitrary (rescaled) 
density and current profiles p,j. Therefore to observe pOOp an average current j over a long time t one should have 



F{j) = lim 



p{x,T),j(x.T) Jo 



dr dx- 



j{x,T) +D{p{x,t)) 



dp{x 



dx 



2a{p{x,T)) 



with the constraint that 



jL = — / j{x,T)dT 



(158) 



(159) 



Comparing (|149p and (|158p we see that the two expressions coincide when the optimal p, j in (|158p are independent 
of the time r. Therefore the additivity principle gives the large deviation function F{j) of the current only when the 
optimal profile in (11581) i s time independent. 

Bertini et al [lOll Il02| pointed out that it can happen, for some a{p) and D{p), that the expression (|151|152p of 
the large deviation function F{j) is non-convex (as it should) in which case the expression (|149p is no longer valid 
(and the prediction (|149p becomes [101, 102] simply an upper bound of F{j)). This is because the optimal profile 
in ()158p is no longer cons tant in time. When this optimal profile is time dependent, one has to solve a much harder 
optimisation problem [sl, Il03l | than (I149P . 

Restrictions on < j(p ) and D{p) for p49p to be valid have been given in [l02j | and there are cases , such as the 
WASEP on a ring [So . Il03| | , for which by varying A or the asymmetry one can observe a phase transition between a 
phase where the optimal profile in p58p is constant in time and a phase where it becomes time dependent. 
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XIV. THE MATRIX APPROACH FOR THE ASYMMETRIC EXCLUSION PROCESS 



The matrix ansatz of section |Vl] (which gives the weights of the microsc opic configuratioiis in the st e ady state) has 
been generahzed to describe the steady state of several other systems [H ES, [lol [loi, [loi, [lol [Tol EM, [HO, [HI 
[m, El, im im im ins ini IHil, [in [m [m , with of course modified algebraic rules for the matrices 
the vectors {W\ and \V). 

For example for the asymmetric exclusion process (ASEP), for which the definition is the same as the SSEP except 
that particles jurnp at r ate 1 to their right and at rate g 7^ 1 to their left it the target site is empty (see figure |4|) , one 
can show IllSj that the weights are still given by (|48p with the algebra (jSOp replaced by 

DE - qED = D + E (160) 
{W\{aE -jD) = {W\ (161) 
{(3D^SE)\V) = \V) . (162) 

One should notice that for the ASEP, the direct approach of calculating, as in section |Vl the steady state properties 
by writing the time evolution leads nowhere. Indeed (I40p becomes 

d{ri) 
dt 

d{n) 

dt 

d{TL) 

dt 



=q; - (a + 7 + 1)(ti) + q{T2) + (1 - q){TiT2) 

- (1 + q){n) + q{n+{) - (1 - q){{n^in) - {nn+i)) (163) 

= -{q + P + &){TL)+5-{\- q){TL-iTL) 



and the equations which determine the one-point functions are no longer closed. Therefore all the correlation functions 
have to be determined at the same time and this is what the matrix ansatz does. Alternative combinatorial 
me thods to c alculate the steady state weights of exclusion processes with open boundary conditions have been obtained 
in [l25l . fT26t . 

The large deviation function T of the density defined by ([67]) has been calculated for the ASEP p?! . Is^ . [53l | by an 
extension of the approach of sections IVIII and IVIIII (see section IXVip . 



XV. THE PHASE DIAGRAM OF THE TOTALLY ASYMMETRIC EXCLUSION PROCESS 

The last three sections IXVIIXVlTl present, as examples, three results which can be obtained rather easily for the 
totally asymmetric case (TASEP) i.e. for g = (in the particular case where particles are injected only at the left 
boundary and removed only at the right boundary i.e. when the input rates 7 = (5 = 0). In this case the algebra 
(|160p becomes 

DE = D + E (164) 
{W\aE={W\ (165) 
PD\V) = \V) (166) 

As for the SSEP the average current (J) is still given in terms of the vectors {W\, V) and of the matrices D and E by 

_ {W\{D + E)^-'\V) 

{W\(D + E)L\V) ■ ^ ' 

However as the algebraic rules have changed, the expression of the current is different for the SSEP and the ASEP. 
From the relation DE = D + E it is easy to prove by recurrence that 

DFiE)=F{l)D + E^^^^l^I^ 

E — 1 

for any polynomial F{E) and 



25 



Using the fact that 

(W|£;™i?"|y) 1 1 



one gi 



ets 



{W\V) a" /?" ' 



1 1 

{W\{D + E)^\V) _ ^ p{2N~l~p)\ ^^gg^ 



{W\V) m{N-p)l 1_1 

a 13 



For large N this sum is dominated either by p ~ 1, or p '--^ A'^ depending one the values of a and [3 and one obtains 

{W\{D + E)^\V) 



4^ if a > i and /3 > i 

2 2 



1 - /3)]-^ if < a and /? < i (169) 



[q;(1 — a)] ^ if /? > a and a < - 



This leads to three different expressions of the current (|167p for large L corresponding to the three different phases: 

• the low density phase {(3 > a and a < ^) where (J) = a{l — a) 

• the high density phase {a> (3 and /? < ^) where (J) = — (3) 

• the maximal current phase (a > | and /? > i) where (J) = -j 

which is the exact phase diagram of the TASEP [§, [l^, [ll|, Uli • The existence of phase transitions [H, [H, [63, [sl] in 
these driven lattice gases is one of the striking properties of non-equilibrium steady states, as it is well known that 
one dimensional systems at equilibrium with short range interactions cannot exhibit phase transitions. 

XVI. ADDITIVITY AND LARGE DEVIATION FUNCTION FOR THE TASEP 

Let us now see how the additivity relation ([TO]) can be generalized for the TASEP in order to obtain the large 
deviation functional of the density. For the algebra (|164m66|) . if one inroduces the following eigenvectors 

{p\E=-{p\ ; D\p) = -^\p) (170) 

P 1 - P 

it is clear that 

{W\ = {pa\ ; \V)^\p,) (171) 

with Pa = a and pb = I — p. Note that in general {pa\pb) ^ even when pa ^ pb- Now one can prove, as in (|60p. that 
for Pb < Pa 

{Pa\YlY2\pb) _ / dp [Pa-Pb) {pa^P) {p\Y2\pb) ^^^2) 



{Pa\pb) J 2iTT (pa - p)ip - Pb) {Pa\p) {p\Pb) 

Pb<\p\<pa 

where Yi and Y2 are arbitrary polynomials in D and E. 

Proof of (|172p : Any polynomial Y of the operators D and E can be written, using DE = D + E, as 

Y = Y,an,n'E"D"' (173) 

n,n' 

by pushing all the D's to the right and all the i?'s to the left. Therefore to prove (|172p it is sufficient to do it for Yi 
and Y2 of the form 

Yi = E'^'D< , Y2 ^ E''^ D"^'^ . (174) 
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If n'l = or 712 = 0, the identity (|172p is easy to check. Then one can prove it by recursion on n'l +712: if Yi = ZiD 
and Y2 = EZ2, and one assumes that the identity ()172p is vahd for Z1DZ2 and Z1EZ2, the l.h.s. of (|172p can be 
written as (|164|) 



{pa\ZiDEZ2\pb) {pa\Zi{D + E)Z2\pb) 



{Pa\Pb) 



{Pa\pb) 



dp {Pa - Pb) 



2i7r {pa - p){p- Pb) 



1 1 



P 1 - P, 



{Pa\Zl\p) {p\Z2\pb) 
{Pa\p) {P\pb) 



Pb<\p\<Pa 

and simply beacuse 1/p + 1/(1 — p) = l/{p{l — p)) this becomes 

dp {.Pa - Pb) 1^ {Pa\Zl\p) {p\Z2\pb) ^ I ^ {pa~Pb) {pa\ZiD\p) {p\EZ2\pb} 

2ilT {pa - p){p - Pb) p{l - p) {Pa\p) {p\Pb) J 2iTr {pa - p){p - Pb) {Pa\p) {p\Pb) 

Pb<\p\<Pa Pb<\p\<Pa 



which is the r.h.s. of p72p . 

We are now going to see, as an example, how the large deviation function T of the density can be derived 
for the TASEP from (|172p when pa > pb- If one defines K{pa, pb) by 



^ r li {Pa\{D + E)^\pb) 
K{pa,Pb) = lim -log — — 

L^oo L [PalPb) 



(175) 



one can easily check from (|169p that for pb{= 1 — P) < Pa{— a) 



K{pa,Pb)^- max log(p(l-p)) 

Pb<P<Pa 



(176) 



Using a saddle point method in (|172p when Yi and I2 are sums of long products of D's and E's, one gets 

{Pa\YlY2\pb) . {Pa\Y,\F) {F\Y2\pb) 

~ mm 

{Pa\pb) Pb<F<p^ {pa\F) (F\pb) 



(177) 



(Note that in applying the saddle point method, one needs to find the maximum F over the circular integration 
contour. This maximum is at the same time a minimum when F varies along the real axis). Then as for a system of 
large size L (|175ll76p one has 



{Pa\Y\pb) 



{pa\{D + E)L\pb) 



{Pa\Y\pb) 
{Pa\Pb) 



(178) 



One can of course repeat p77p several times to relate a large system of size L to its subsystems (as long as these 
subsystems are large). Therefore one gets 



{Pa\YiY2...Yk\pb) _ g-K(p„p,) L ^-^ TT 

{pa\{D + L)^\pb) p^>F,>F2...>F„_,>p, 



{Fi^l\Yi\Fi) K(Fi-i.F,) I 

{F,^i\{D + Ey\F,) 



(179) 



where Fq = Pa, Fn = pb and I — L/n. If Yi is the sum of the matrix elements of all configurations of / sites with Ipi 
particles one gets for the large deviation function !Fn defined in 



J'7iipi,P2,--Pn\Pa,Pb) ^ K{pa,Pb) + ^ ^^ax 

Pb=Fo<Fi..<Fi<..<F„=Pa n 



1 " 

-y^J^i{p,\F,.i,F,) - K{F,^i,F,) . 



(180) 



For large n, almost all the differences i^i-i — Fi are small, so that 



Tiip,\F,^i,F,) ~ T,ip,\F,,F,) = log § + (1 - P^) log ^— ^ = B(p„ F,) 



since when the two densities Fi-i and Fi are equal, the steady state measure is Bernoulli and this leads to 



H{pi^)}) 



max 

Pb<r<Pa 



log (r(l - r)) 



sup / dx 
F Jo 



B 



{p{x),F{x))+\og(^F{x){l-F{x)))] , (181) 
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which is the expression of the large deviation function of the density of the TASEP (and also of the ASEP [sl, ISSj ') 

for Pa> Pb- 

For Pa > Pb, one can also obtain [23, [52, \EM ^^is large deviation function, starting from a relation similar to (|172p 
obtained by deforming the circular contour to insure that pb remains inside and pa outside. One can note that when 
Yi and Y2 are polynomials in D and E, all the matrix elements in (|172[) are rational functions of pa and ph which can 
be easily anaytically continued from the case pa > Pb to the case pa < Pb- The result is [sH [53l| 



J'{{p{x)})= inf 

0<y<l 



\x (^Bip{x),pa) + \ogP^^L_PA 

Pb(l - Pb) 

(J) 



dx B{p{x),pb) + log 



(182) 



For the TASEP one knows [l^, [Tl| that along the line Pa = 1 — Pfc < 5 there is a first order phase tr ansition line. 
Along this line (J) = pa(l — pa) = pb{^ — Pb) and the typical configurations Pz{x) are schocks |ll3l . IllSl [l2^7j located 
at arbitrary positions z beween a region of density pa and a region of density 1 — pa- 



Pzix) 



Pa for < a; < z 
1 — Pa for z < X < 1 



(183) 



For all these profiles Pzix), the functional (|182|) vanishes. It is also easy to check that !F{p{x)) > for a profile of the 
form 



p{x) = apz{x) + (1 - a)pz'{x) 



(184) 



and this shows that J- is non-convex. Therefore in contrast to equilibrium systems, the functional !F{p{x)) may be 
non-convex in non-equilibrium steady states. 



XVII. CORRELATION FUNCTIONS IN THE TASEP AND BROWNIAN EXCURSIONS 



In this last example, we will see that the fiuctuations of the density are non-Gaussian in the maximal current phase 
of the TASEP. In this maximal current phase (a > ^ and (3 > ^) one can show [s^l, using the matrix ansatz, that the 
correlation function of the occupations of k sites at positions ii = Lxi, 1-2 = Lx2, ■■■ik = Lxk with xi < X2 < ... < Xk 
are given by 

( ("^- ~ ^) • ■ ■ ("^- - ^) ) ^ ¥lkdx,'^..dxM'''^ ■ ■ ■ ^^"'■^) ' ^'^^^ 

where y{x) is a Brownian excursion between and 1 (a Brownian excursion is a Brownian path constrained to 
y{x) > for < a; < 1 with the boundaries y(0) = y(l))- The probability P(yi . . .yk;xi . . .Xk) of finding the 
Brownian excursion at positions yi ... yk for < xi <...< Xk < I is 

N h^iiyi) 9x2-xiiyi,y2) ■■ -gx^-Xk-iivk-i^yk) hi^^^{yk) 

P[yi . . .yk]Xi . . .Xk) = -j= , 



where hx and gx are defined by 



2y .-y'lx 



gx{y.y') - [e-^v-vr/^^ _ ^-(v+v'?/^'^ . 

One can derive easily (|185p in the particular case a — (3 = \ using a representation of (I164p which consists of two 
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infinite dimensional bidiagonal matrices 



Z? = ^|n)(n| + |n)(,^+l| = 



ii;-^|n)(n+l| + |n)(n| = 



n>l 



/I 1 

110 






110 
1 1 



V 

/I 

1 1 

1 1 

1 













with 

(W^| = (l| = (1,0,0...) 
(F| = (1| = (1,0,0...). 

With this representation one can write (VFKI? + E)^\V) as a sum over a set Ai^ of one dimensional random walks w 
of L steps which remain positive. Each walk w is defined by a sequence (ni(w)) of L — 1 heights {ni{w) > 1) (with 
at the boundaries no{w) — nL{w) — 1 and the constraint — n^j < 1) : 

{W\{D + E)^\V) ^ n{w), 



where 



^('"^) — JJ^ ''^('^i-ii '^i) with v(n,n'^ 



2 if \n-n'\ =0 
1 if I n — n' I = 1 . 



One has v[n, n') — {n\D + E\n') since D + E has a tridiagonal form 

/ 2 1 (0)\ 
1 '•. '•• 



D + E 



1 



V(0) 1 2 J 



Then from the matrix expression one gets (t;) et {riTj 



where v(w) is the probability of the walk w induced by the weights O : 



v{w) 



(186) 



The expression (|186p is the discrete version of (llSSp . The result (|185|) can be extended [5j| to arbitrary values of a 
and p in the maximal current phase (i.e. for a > 1/2 and /3 > 1/2). 

From this link between the density fluctuations and Brownian excursions, one can show that, for a TASEP of L 
sites, the number N of particles beween sites Lxi and Lx2, has non-Gaussian fluctuations in the maximal current 
phase: if one defines the reduced density 



N~L{x2-xi)/2 



(187) 
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one can show ^53| that for large L 



P{fi) = dyi I dy2 = 
Jo Jo ^J2-K{xi~xx) 

This contrasts with the Gaussian fluctuations of the density (|24p at equihbrium. According to numerical simulations 
[5^ the distribution (|187ll88p (properly rescaled) of the fluctuations of the density remains valid for more general 
driven systems in their maximal current phase. Of course proving it in a more general case is an interesting open 
question. 



X2 — Xl 



(188) 



XVIII. CONCLUSION 



These lectures have presented a number of recent results concerning the fluctuations and the large deviation 
functions of the density or of the current in non-equilibrium steady states. 

For general diffusive systems, the macroscopic fluctuation theory [4l|, [H, IH, [H, llOll . Il02l | discussed in sec- 
tion [IX] gives a framework to calculate the large deviation of the density T{p{x)) leading to equations (I96I95I94P 
which are in general difficult to solve. One can however check that the expressions (|71l72l78p obtained [2^, 
sections IVHIVIIII by the matrix ansatz do solve these equations. So far the large deviation functional is only known 
for very few examples [2^, [2^, [13, There remains a lot to be done to understand the general properties of 

the functional !F{p{x)). For example, with Thierry Bodineau we tried, so far without success, to calculate J-{p{x)) 
(|96l95l94p for general D and a in powers of Pa — Pb- Also for the SSEP, we did not succeed to obtain the large 
deviation functional of the density T{p{x)) for the A— measure defined in (|129p . Other situations which wo uld 
be interesting to conside r are the cases of several species of particles, several conserved quantities [l28l Il29l | or 
non-conserved quantities [l30j | . 

For driven diffusive systems, the situation is worse: so far T{p{x)) is only known for the ASEP [s^, 113 smii, 
to my knowledge, there does not exist so far a general theory to calculate this large deviation functional for general 
driven diffusive systems |l3l| . In contrast to equilibrium systems, for the ASEP, the large deviation functional may 
be non-convex fsection [XVip or fluctuations of the density may be non-Gaussian (section [XVIip . 



For current fluctuations in diffusive systems, the additivity principle gives explicit expressions (section IXIIip 
of the large deviation function of the current for general diffusive systems. In some cases, however, these expressions 
are not valid, when the profile to produce a deviation of current becomes time dependent. In such cases the calculation 
of the large deviation function F{j) of the current is much harder [s^. So far the predictions (|149ll51ll52ll54j|156p of 
the additivity principle remain to be checked in specific examples: even for the SSEP, only the first four cumulants are 
known to agree with (|156p but a direct calculation of F{j) for the SSEP is, to my knowledge, still missing. It wo uld 
be nice to see whether this could be done by a Bethe ansatz calculation for the SSEP with open boundaries |l32l | . It 
would also be useful to test the predictions of the additivity principle on other diffusive systems and to try to ex - 
tend them to more complicated situations, in particular when there are more than a single conserved quanity |l28l . ll2i^ . 

Concerning the fluctuations or the l arge deviations of the curr ent of driven diffusive systems, there has been 
lots of progress over the last ten years [lllIlMIllEIIlEIIllIlll. On the infinite line exact resuhs for the TASEP 
and the PNG (polynuclear growth model) lead to a universal distribution of c urrent characteristic of the KPZ 
universality class. On the ring too, Bethe ansatz calculations [65, 95, 96, 13l [Tiol [Till . [lil [Tisl . [Til [Tisl [l4^ [T47t. 
allow to calculate the large deviation function of the current which exhibits a universal shape in the scaling regime. 
For driven diffusive systems, however, there is not yet a general approach allowing to calculate the large deviation 
func tion or the fluctuations of the current for all geometries, including finite systems with open boundary conditions 
[l32| . Of course it would be nice to extend the macroscopic fluctuation theory to get a framework allowing to 
calculate both the large deviation functions of the current and of the density for general driven diffusive systems. 
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years. I would like to thank C. Appert, T. Bodineau, M. Clincy, B. Dougot, C. Enaud, M.R. Evans, S. Goldstein, V. 
Hakim, S.A. Janowsky, C. Landim, J.L. Lebowitz, K. Mallick, D. Mukamel, S. OUa, V. Pasquier, P.E. Roche, E.R. 
Speer with whom I had the pleasure and the privilege to work on these topics. 



[1] F. Bonetto, J.L. Lebowitz, L. Rey-Bellet, Fourier's law: a challenge to theorists, Mathematical Physics 2000, 128-150, 

Imperial College Press (2000). math-ph/0002052 
[2] S. Lepri, R. Livi, A. Politi, Thermal conduction in classical low- dimensional lattices, Phys. Rep. 377, 1-80 (2003) 
[3] J. P. Eckmann, C.A. Pillet, L. Rey-Bellet, Entropy production in nonlinear, thermally driven Hamiltonian systems J. Stat. 

Phys. 95, 305-331 (1999) 
[4] S.L Sasa , H. Tasaki, Steady state thermodynamics J. Stat. Phys. 125 125-227 (2006) 

[5] D. Ruelle, Conversations on nonequilibrium physics with an extraterrestrial, Physics Today 57, No 5, 48-53 (2004) 
[6] D. Ruelle, Smooth dynamics and new theoretical ideas in nonequilibrium statistical mechanics, J. Statist. Phys. 95, 393-468 
(1999) 

[7] S. Katz, J. Lebowitz, H.Spohn, N on- equilibrium steady states of stochastic lattice gas models of fast ionic conductors, J. 

Stat. Phys. 34, 497-537 (1984) 
[8] E.D. Andjel, Invariant-measures for the zero range process Ann. Prob. 10 525-547 (1982) 

[9] B. Derrida, E. Domany, D. Mukamel, An exact solution of a one- dimensional asymmetric exclusion model with open 

boundaries, J. Stat. Phys. 69, 667-687 (1992) 
[10] B. Derrida, M.R. Evans, V. Hakim, V. Pasquier, Exact solution of a Id asymmetric exclusion model using a matrix 

formulation, J. Phys. A 26, 1493-1517 (1993) 
[11] G. Schiitz, E. Domany, Phase transitions in an exactly soluble one- dimensional asymmetric exclusion model, J. Stat. Phys. 

72, 277-296 (1993) 

[12] B. Derrida, S.A. Janowsky, J.L. Lebowitz, E.R. Speer, Exact solution of the totally asymmetric simple exclusion process 

- shock profiles J. Stat. Phys. 73 813-842 (1993) 
[13] M.R.Evans, S.N. Majumdar, R.K.P. Zia, Factorized steady states in mass transport models J. Phys A 37 L275-L280 (2004) 
[14] J. P. Eckmann, L.S. Young, Temperature profiles in Hamiltonian heat conduction Europhys. Lett. 68 790-796 (2004) 
[15] J. P. Eckmann, L.S. Young, Nonequilibrium energy profiles for a class of 1-D models Comm. Math. Phys. 262 237-267 

(2006) 

[16] F. Bonetto, J.L. Lebowitz, J. Lukkarinen, Fourier's law for a harmonic crystal with self-consistent stochastic reservoirs 

J. Stat. Phys. 116 : 783-813 (2004) 
[17] Y. Srebro, D. Levine, Exactly solvable model for driven dissipative systems Phys. Rev. Lett. 93 240601 (2004) 
[18] O. Narayan, S. Ramaswamy, Anomalous heat conduction in one- dimensional momentum- conserving systems Phys. Rev. 

Lett. 89 200601 (2002) 

[19] S. Lepri, R. Livi, A. Politi, On the anomalous thermal conductivity of one- dimensional lattices Europhys. Lett. 43 271-276 
(1998) 

[20] G. Eyink, J.L. Lebowitz, H. Spohn, Lattice gas models in contact with stochastic reservoirs: local equilibrium and relaxation 

to the steady state. Comm. Math. Phys. 140 119-131 (1991) 
[21] C. Kipnis, S. OUa, S. Varadhan, Hydrodynamics and large deviations for simple exclusion processes, Commun. Pure Appl. 

Math. 42, 115-137 (1989) 
[22] H. Spohn, Large scale dynamics of interacting particles ( Springer- Verlag, Berlin, 1991) 

[23] T. Liggett, Stochastic interacting systems: contact, voter and exclusion processes. Fundamental Principles of Mathematical 

Sciences, 324 Springer- Verlag, Berlin, (1999) 
[24] C. Kipnis, C. Landim, Scaling limits of interacting particle systems. Springer (1999) 

[25] B. Derrida, J. L. Lebowitz, E. R. Speer, Free energy functional for nonequilibrium systems: an exactly solvable case, Phys. 
Rev. Lett. 87 150601 (2001) 

[26] B. Derrida, J.L. Lebowitz, E.R. Speer, Large deviation of the density profile in the steady state of the open symmetric 

simple exclusion process, J. Stat. Phys. 107 599-634 (2002) 
[27] C. Enaud, B. Derrida, Large deviation functional of the weakly asymmetric exclusion process, J. Stat. Phys. 114 537-562 

(2004) 

[28] E.D. Andjel , M. D. Bramson, T.M. Liggett, Shocks in the asymmetric exclusion process, Prob. Th. and Rel. Fields 78 
231-247 (1988) 

[29] P.A. Ferrari, C. Kipnis, E. Saada, Microscopic structure of traveling waves in the asymmetric simple exclusion process, 

Annals of Probability 19 226-244 (1991) 
[30] P.A. Ferrari, Shock fluctuations in asymmetric simple exclusion, Prob. Th. and Rel. Fields 91 81-101 (1992) 
[31] J. Krug, Boundary-induced phase-transitions in driven diffusive systems, Phys. Rev. Lett. 67 1882-1885 (1991) 
[32] R. Bundschuh Asymmetric exclusion process and extremal statistics of random sequences PHys. Rev. E 65 031911 (2002) 
[33] L.B. Shaw, R.K.P. Zia, K.H. Lee Totally asymmetric exclusion process with extended objects: A model for protein synthesis 

Phys. Rev. 68 021910 (2003) 

[34] P.M. Richards, Theory of one- dimensional hopping conductivity and diffusion Phys. Rev. B 16 1393-1409 (1977) 



31 



D. Chowdhury, L. Santen, A. Schadschneider, Statistical physics of vehicular traffic and some related systems, Phys. Rep. 
329 199-329 (2000) 

T. Halpin-Hcaly, Y.C. Zhang, Kinetic roughening phenomena, stochastic growth directed polymers and all that Phys. Rep. 
254 215-415 (1995) 

B. Widom, J.L. Viovy, A.D. Defontaines, Repton model of gel- electrophoresis and diffusion J. de Physique I 1 1759-1784 
(1991) 

M. Kardar, G. Parisi, Y.C. Zhang, Dynamic scaling of growing interfaces Phys. Rev. Lett. 56 889-892 (1986) 

J.E. Santos and G. M. Schiitz, Exact time- dependent correlation functions for the symmetric exclusion process with open 

boundary 64, 036107 (2001) 

S.F. Edwards, D.R. Wilkinson, The surface statistics of a granular aggregate Proc. Roy. Soc. A 381 17-31 (198a)2 

L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, C. Landim, Fluctuations in stationary non equilibrium states of 

irreversible processes, Phys. Rev. Lett. 87 040601 (2001) 

L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, C. Landim, Macroscopic fluctuation theory for stationary non 
equilibrium states, J. Stat, Phys. 107, 635-675 (2002) 

L. Bertini, A. De Solo, D. Gabrielli, G. Jona-Lasinio, C. Landim, Large deviations for the boundary driven symmetric 
simple exclusion process. Math. Phys. Analysis and Geometry 6, 231-267 (2003) 

L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, C. Landim, Minimum dissipation principle in stationary non 
equilibnum states, J. Stat. Phys. 116 831-841 (2004) 

D.J. Evans, E.G.D. Cohen, G.P. Morriss, Probability of second law violations in shearing steady states, Phys. Rev. Lett. 
71, 2401 (1993) 

G. Gallavotti, E.D.G. Cohen, Dynamical ensembles in stationary states, J. Stat. Phys. 80, 931-970 (1995) 
D.J. Evans, D.J. Searles, The Fluctuation Theorem, Advances in Physics 51, 1529-1585 (2002) 
J. Kurchan, Fluctuation Theorem for stochastic dynamics, J. Phys. A31 3719, (1998) 

J.L. Lebowitz, H. Spohn, A Gallavotti-Cohen type symmetry in the large deviation functional for stochastic dynamics, J. 
Stat. Phys. 95, 333-365 (1999) 

B. Derrida, B. Dougot, P.-E. Roche, Current fluctuations in the one- dimensional symmetric exclusion process with open 
boundaries, J. Stat. Phys. 115, 717-748 (2004) 

T. Bodineau, B. Derrida, Current fluctuations in nonequilibriurn diffusive systems: An additivity principle Phys. Rev. 
Lett. 92, 180601 (2004) 

B. Derrida, J.L. Lebowitz, E.R. Speer, Exact free energy functional for a driven diffusive open stationary nonequilibrium 
system, Phys. Rev. Lett. 89 030601 (2002) 

B. Derrida, J.L. Lebowitz, E.R. Spcer, Exact large deviation functional of a stationary open driven diffusive system: the 
asymmetric exclusion process, J. Stat. Phys. 110 775-810 (2003) 

B. Derrida, C. Enaud, J.L. Lebowitz, The asymmetric exclusion process and Brownian excursions, J. Stat. Phys. 115 
365-382 (2004) 

J. P. Eckmann, C.A. Pillet, L. Rey-Bellet, Non-equilibrium statistical mechanics of anharmonic chains coupled to two heat 
baths at different temperatures. Comm. Math. Phys. 201, no. 3, 657-697 (1999) 

T. Bodineau, B. Derrida, Cumulants and large deviations of the current through non- equilibrium steady states, Submited 
to CRAS 2007 

C. Macs, F. Rcdig F, A.V. Moffaert On the definition of entropy production, via examples J. Math. Phys. 41 1528-1554 
(2000) 

U. Seifert, Entropy production along a stochastic trajectory and an integral fluctuation theorem Phys. Rev. Lett. 95 040602 
(2005) 

R. Ellis, Entropy, large deviations, and statistical mechanics. Reprint of the 1985 original. Classics in Mathematics 

Springer- Vcrlag, Berlin, (2006) 

A. Dcrnbo, O. Zcitouni, Large deviations techniques and applications. Second edition Applications of Mathematics, 38. 
Springer- Ver lag, (1998) 

M.R. Evans, D.P. Foster, C. Godreche, D. Mukamel, Spontaneous symmetry-breaking in a one- dimensional driven diffusive 
system Phys. Rev. Lett. 74 208-211 (1995) 

K. Mallick, Shocks in the asymmetry exclusion model with an impurity J. Phys. A 29 5375-5386 (1996) 

M.R. Evans, Y. Kafri, H.M. Koduvely, D. Mukamel, Phase separation in one- dimensional driven diffusive systems Phys. 

Rev. Lett. 80 425-429 (1998) 

A. B. Kolonicisky, G.M. Schutz, E.B.Kolomeisky, J. P. Straley, Phase diagram of one- dimensional driven lattice gases with 
open boundaries J. Phys. A 31 6911-6919 (1998) 

B. Derrida, An exactly soluble non-equilibrium system: the asymmetric exclusion process Phys. Rep. 301 , 65-83 (1998) 
V. Popkov, G.M. Schutz Steady-state selection in driven diffusive systems with open boundaries Europhys. Lett. 48 257-263 

(1999) 

M.R. Evans, Phase transitions in one- dimensional nonequilibrium systems, Braz. J. Phys. 30, 42-57 (2000) 

Y. Kafri, E. Lcvinc, D. Mukanicl, G.M. Schiitz, J. Torok, Criterion for phase separation in one- dimensional driven 

systems, Phys. Rev. Lett. 89, 035702 (2002) 

A. Parmeggiani, T. Franosch, E. Frey, Phase coexistence in driven one- dimensional transport Phys. Rev. Lett. 90 086601 

(2003) 

M. Chncy, B. Derrida, M.R. Evans, Phase transitions in the ABC model, Phys. Rev. E 67, 066115 (2003) 

G. Fayolle, C. Furtlehner, Dynamical windings of random walks and exclusion models. Part /: Thermodynamic limit in 



32 



Z(2) J. Stat. Phys. 114 229-260 (2004) 
[72] H. Spohn, Long range correlations for stochastic lattice gases in a nori- equilibrium steady state, J. Phys. A 16 4275-4291 
(1983) 

[73] R. Schmitz, E.G.D. Cohen, Fluctuations in a fluid under a stationary heat-flux .1. General theory J. Stat. Phys. 39 285-316 
(1985) 

[74] J.R. Dorfman, T.R. Kirkpatrick, J.V. Sengers, Generic long-range correlations in molecular fluids Annual Review of 

Physical Chemistry 45 213-239 (1994) 
[75] J.M. O. de Zarate, J.V. Sengers, On the physical origin of long-ranged fluctuations in fluids in thermal nonequilihrium 

states J. Stat. Phys. 115 1341-1359 (2004) 
[76] B. Derrida, C. Enaud, C. Landirn, S. OUa, "Fluctuations in the weakly asymmetric exclusion process with open boundary 

conditions" J. Stat. Phys. 118, 795-811 (2005) 
[77] B. Derrida, J.L. Lebowitz, E.R. Speer, Entropy of open lattice systems, to appear in J. Stat. Phys. 2007 
[78] V. Hakim, J. P. Nadal, Exact results for 2d directed animals on a strip of finite width J. Phys. A 16, L213-L218 (1983) 
[79] I. Affleck, T. Kennedy, E.H. Lieb, H. Tasaki, Valence bond ground-states in isotropic quantum antiferromagnets Comm. 

Math. Phys. 115 477-528 (1988) 
[80] A. Klumper, A. Schadschneider, J. Zittartz, Equivalence and solution of anisotropic spin-1 models and generalized T-J 

Fermion models in one dimension J. Phys A 24 L955-L959 (1991) 
[81] S. Pilgram, A.N. Jordan, E.V. Sukhorukov, M. Buttiker, Stochastic path integral formulation of full counting statistics, 

Phys. Rev. Lett. 90, 206801 (2003) 
[82] A.N. Jordan, E.V. Sukhorukov, S. Pilgram, Fluctuation statistics in networks: A stochastic path integral approach, J. 

Math. Phys. 45 4386-4417 (2004) 
[83] C. Kipnis, C. Marchioro, E. Presutti, Heat-flow in an exactly solvable model, J. Stat. Phys. 27 65-74 (1982) 
[84] L. Bertini, D. Gabrielli, J. Lebowitz, Large deviation for a stochastic model of heat flow, J. Stat. Phys. 121, no. 5-6, 

843-885 (2005) 

[85] A. Demasi, E. Presutti, E. Scacciatelli, The weakly asymmetric simple exclusion process, Ann. Institut Henri Poincare- 
Prob. Stat. 25, 1-38 (1989) 

[86] R. J. Harris, A. Rakos, G. M. Schiitz, Current fluctuations in the zero-range process with open boundaries, J. Stat. Mech. 
P08003 (2005) 

[87] G. Gallavotti, Chaotic hypothesis: Onsager reciprocity and fluctuation-dissipation theorem, J. Stat. Phys. 84, 899-926, 
(1996) 

[88] G.E. Crooks, Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences 

Phys. Rev. E 60 2721-2726 (1999) 
[89] C. Macs, The fluctuation theorem as a Gibbs property, J. Stat. Phys. 95, 367-392 (1999) 

[90] G. Gallavotti, Entropy production in nonequilibrium thermodynamics: a point of view. Chaos 14, 680-690, (2004) 
[91] P. Gaspard, Fluctuation theorem for nonequilibrium reactions J. Chem. Phys. 120 8898-8905 (2004) 
[92] J. Farago, Power fluctuations in stochastic models of dissipative systems, Physica A 331, 69-89 (2004) 
[93] R. J. Harris, A. Rakos, G. M. Schiitz, Breakdown of Gallavotti-Cohen symmetry for stochastic dynamics, Europhysics 
Letters, 75, 227-233 (2006) 

[94] P. Visco, Work fluctuations for a Brownian particle between two thermostats, J. Stat. Moch. P06006 (2006) 
[95] B. Derrida, J.L. Lebowitz, Exact large deviation function in the asymmetric exclusion process, Phys. Rev. Lett. 80, 
209-213 (1998) 

[96] B. Derrida, C. Appert, Universal large deviation function of the Kardar-Parisi- Zhang equation in one dimension J. Stat. 
Phys. 94, 1-30 (1999) 

[97] C. Giardina, J.Kurchan, L. Politi, Direct evaluation of large-deviation functions Phys. Rev. Lett. 96 120603 (2006) 
[98] H. Lcc, L. S. Lcvitov, A. Yu. Yakovots, Universal statistics of transport in disordered conductors, Phys. Rev. B 51, 
4079-4083 (1995) 

[99] Y.M. Blanter, M. Buttiker, Shot noise in mesoscopic conductors, Phys. Rep. 336, 1-166 (2000) 
[100] F. van Wijland, Z. Racz, Large deviations in weakly interacting boundary driven lattice gases, J. Stat. Phys. 118, 27-54 
(2005) 

[101] L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, C. Landim, Current fluctuations in stochastic lattice gases, Phys. 
Rev. Lett. 94 030601 (2005) 

[102] L. Bertini, A. Dc Sole, D. Gabrielli, G. Jona-Lasinio, C. Landirn, Non equilibrium current fluctuations in stochastic lattice 

gases, J. Stat. Phys. 123, no. 2, 237-276 (2006) 
[103] T. Bodineau, B. Derrida, Distribution of current in nonequilibrium diffusive systems and phase transitions, Phys. Rev. E 

72, no. 6, 066110 (2005) 

[104] B. Derrida, M.R. Evans, D. Mukamel, Exact diffusion constant for one dimensional asymmetric exclusion models J. Phys. 
A26 4911-4918 (1993) 

[105] S. Sandow, Partially asymmetric exclusion process with open boundaries, Phys. Rev. E 50 2660-2667 (1994) 

[106] B. Derrida, M.R. Evans, K. Mallick, Exact diffusion constant of a one dimensional asymmetric exclusion model with open 

boundaries J. Stat. Phys. 79 833-874 (1995) 
[107] R.B. Stinchcombe, G.M. Schutz, Application of operator-algebras to stochastic dynamics and the Heisenberg chain Phys. 

Rev. Lett. 75 140-143 (1995) 

[108] B. Derrida, K. Mallick, Exact diffusion constant for the one- dimensional partially asymmetric exclusion model J. Phys. 
A Phys. 30 1031-1046 (1996) 



33 



[109] F.H.L. Essler, V. Rittenberg, Representations of the quadratic algebra and partially asymmetric diffusion with open 

boundaries , J. Phys. A 29 3375-3407 (1996) 
[110] H. Hinrichsen, S. Sandow, I. Peschel, On matrix product ground states for reaction-diffusion models, J. Phys. A 29 

2643-2649 (1996) 

[111] K. Mallick, S. Sandow, Finite- dimensional representations of the quadratic algebra: Applications to the exclusion process, 

J. Phys. A 30 4513-4526 (1997) 
[112] E.R. Speer, Ftnite-dimensional representations of a shock algebra, J. Stat. Phys. 89, 169-175 (1997) 

[113] B. Derrida, J.L. Lebowitz, E.R. Speer, Shock profiles for the partially asymmetric simple exclusion process J. Stat. Phys. 
89, 135-167 (1997) 

[114] K. Krebs, S. Sandow, Matrix product eigenstates for one- dimensional stochastic models and quantum spin chains, J. Phys. 
A 30 3165-3173 (1997) 

[115] B. Derrida, S. Goldstein, J.L. Lebowitz, E.R. Speer, Shiff equivalence of measures and the intrinsic structure of shocks in 
the asymmetric simple exclusion process J. Stat. Phys. 93, 547-571 (1998) 

[116] T. Sasamoto, One-dimensional partially asymmetric simple exclusion process with open boundaries: orthogonal polynomi- 
als approach, J. Phys. A 32 7109-7131 (1999) 

[117] F.C. Alcaraz , S. Dasmahapatra, V. Rittenberg, N-species stochastic models with boundaries and quadratic algebras, J. 
Phys. A 31 845-878 (1998) 

[118] R.A. Blythe, M.R. Evans, F. Colaiori, F.H.L. Essler, Exact solution of a partially asymmetric exclusion model using a 

deformed oscillator algebra, J. Phys. A 33 2313-2332 (2000) 
[119] H. Hinrichsen, N on- equilibrium critical phenomena and phase transitions into absorbing states. Advances in Physics 49 

815-958 (2000) 

[120] A.P. Isaev, P.N. Pyatov, V. Rittenberg, Diffusion algebras, J. Phys. A 34 5815-5834 (2001) 

[121] C. Boutillier, P. Frangois, K. Mallick, S. Mallick, A matrix ansatz for the diffusion of an impurity in the asymmetric 

exclusion process, J. Phys. A 35 9703-9730 (2002) 
[122] F.H. Jafarpour, Matrix product states of three families of one- dimensional interacting particle systems, Physica A 339 

369-384 (2004) 

[123] M. Depken, R. Stinchcombe, Exact Joint Density -Current Probability Function for the Asymmetric Exclusion Process, 

Phys. Rev. Lett. 93, 040602 (2004) 
[124] M. Depken, R. Stinchcombe, Exact probability function for bulk density and current in the asymmetric exclusion process, 

Phys. Rev. E 71, 036120 (2005) 
[125] R. Brak, J.W. Essam Asymmetric exclusion model and weighted lattice paths J. Phys. A 37 4183-4217 (2004) 
[126] E. Duchi, G. SchaefFer, A combinatorial approach to jumping particles J. Comb. Th. A 110 1-29 (2005) 
[127] L. Santen, C. Appert, The asymmetric exclusion process revisited: Fluctuations and dynamics in the domain wall picture 

J. Stat. Phys. 106 187-199 (2002) 
[128] J. Fritz, K. Nagy, S. OUa, Equilibrium fluctuations for a system of harmonic oscillators with conservative noise J. Stat. 

Phys. 122 399-415 (2006) 

[129] G. Basile, C. Bernardin, S. OUa Momentum conserving model with anomalous thermal conductivity in low dimensional 

systems Phys. Rev. Lett. 96 204303 (2006) 
[130] G. Basile, G. Jona-Lasinio, Equilibrium states with macroscopic correlations Int. J. Mod. Phys. B 18 479-485 (2004) 
[131] T. Bodineau, B. Derrida, Current large deviations for asymmetric exclusion processes with open boundaries, J. Stat. Phys. 

123, no. 2, 277-300 (2006) 

[132] J. de Gier, F.H.L. Essler, Bethe ansatz solution of the asymmetric exclusion process with open boundaries Phys. Rev. Lett. 
95 240601 (2005) 

[133] K. Johansson, Shape fluctuations and random matrices Comm. Math. Phys. 209 437-476 (2000) 

[134] M. Prahofer, H. Spohn, Universal distributions for growth processes in l-hl dimensions and random matrices Phys. Rev. 

Lett. 84 240601 (2005) 4882-4885 (2000) 
[135] M. Prahofer, H. Spohn, Scale invariance of the PNG droplet and the Airy process J. Stat. Phys. 108 1071-1106 (2002) 
[136] T. Sasamoto, Spatial correlations of the ID KPZ surface on a flat substrate J. Phys. A 38 L549-L556 (2005) 
[137] H. Spohn Exact solutions for KPZ-type growth processes, random matrices, and equilibrium shapes of crystals Physics A 

369 71-99 (2006) 

[138] A. Borodi, P. L. Ferrari, M. Prahofer, T. Sasamoto, Fluctuation properties of the TASEP with periodic initial configuration 
math-ph/ 0608056 

[139] L.H. Gwa, H. Spohn, Bethe solution for the dynamic-scaling exponent of the noisy Burgers equation Phys. A 46 844-854 
(1992) 

[140] G.M. Schutz, Exact solution of the master equation for the asymmetric exclusion process J. Stat. Phys. 88 427-445 (1997) 
[141] D.S. Lee, D. Kim, Large deviation function of the partially asymmetric exclusion process Phys. Rev. E 59 6476-6482 
(1999) 

[142] V.B. Priezzhev, Exact non stationary probabilities in the asymmetric exclusion process on a ring Phys. Rev. Lett. 91 
050601 (2003) 

[143] J.G. Brankov, V.B. Priezzhev, R.V. Shelest, Ceneralized determinant solution of the discrete-time totally asymmetric 

exclusion process and zero-range process Phys. Rev. E 69 066136 (2004) 
[144] O. Golinelli, K. Mallick, Bethe ansatz calculation of the spectral gap of the asymmetric exclusion process J. Phys. A 37 

3321-3331 (2004) 

[145] A. Rakos, G.M. Schutz, Current distribution and random matrix ensembles for an integrable asymmetric fragmentation 



34 



process J. Sta. Phys. 118 511-530 (2005) 
[146] A.M. Povolotsky , V.B. Priczzhov, Determinant solution for the totally asymmetric exclusion process with parallel update 

J. Stat. Mcch.-Th. Exp. P07002 (2006) 
[147] O. Golinelli, K. Mallick, The asymmetric simple exclusion process: an integrable model for non- equilibrium statistical 

mechanics J. Phys. A 39 12679-12705 (2006) 



35 



